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ABSTRACT 

An  investigation  and  evaluation  is  given  on  the  possible  application 
of  radar-optics  imaging  methods.  The  discussion  concentrates  largely 
on  microwave  holography.  Problem  areas  and  alternatives  are  treated, 
including  a  brief  discussion  on  thinned  array  techniques. 

The  analytical  work  centers  on  array  imaging,  sampled  reconstruc¬ 
tion,  and  sensitivity-noise  effects.  In  addition  to  other  system  considera¬ 
tions,  an  example  is  given  comparing  microwave  and  optical  setups. 
Finally,  an  examination  is  made  of  experimental  back-iq>  in  terms  of 
i*equirements,  advantages,  frequency  selecticm,  and  appnqpriate  optical 
simulation. 
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SECTION  I 

•  BACKGROUND  INFORMATION 

•  HISTORICAL  BACKGROUND 

Studies  were  initiated  under  this  project  very  early  in  1967  to  consider 
radar-optics  methods  using  large  filled  apertures  as  well  as  scattering  center 
map  processing  from  thinned  array  configurations.  Emphasis  was  also  given 
to  wideband  imaging  which  in  particular  corresponded  well  with  the  point  of 
view  adopted  by  the  IDA  subpanel  on  radar  imaging.  This  early  work  on  the 
use  of  radar-optics  analogs  in  image  formation  included  a  preliminary  investi¬ 
gation  of  phase  degradation  due  to  atmosphere  inhomogeneities  on  maintaining 
coherence  across  the  array  in  addition  to  an  establishment  of  the  problem 
^  areas  requiring  study  and  resolution  to  allow  radar-optics  imaging. 

Last  year,  imaging  methods  using  both  thinned  simple  arrays  and 
especially  large  filleo  arrays  were  established  and  suggested  having  potential 
application  to  a  full  range  of  target  sources.  In  this  context  the  possibility  of 
obtaining  direct  real-time  images  using  radar-optics  on  large  arrays  was  to 
be  dealt  with  initially  as  a  feasibility  study.  The  present  report  represents 
the  first  formal  account  of  the  studies  rccomplished  to  date  together  with  an 
indication  of  the  advances  made  along  with  difficulties  and  possibilities 

j 

I  remaining. 

I 

PROBLEM  AREA  AND  IMAGING  ALTERNATIVES 

i 

I  The  target  environment  can  include  both  the  target  source  and  other 

I  background  sources  of  widely  varying  characteristics.  It  is  then  essential 

I  that  high  resolution  be  combined  with  a  real-time  capability  to  produce  images 
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which  can  allow  identification  and  discrimination  in  the  shortest  possible 
time.  Such  a  requirement  dictates  the  use  of  large  filled  arrays. 

Alternatives  using  thinned  array  techniques  or  pure  optical  methods 
could  be  considered. 

Using  relatively  simple  arrays  including  monostatic,  interferometer 
and  trilateration  measurements,  a  number  of  wideband  and  narrowband  in¬ 
verse  methods  have  been  developed  recently  under  tiiis  program  which  create 
maps  of  the  target  as  returns  from  prominent  scattering  centers.  These 
maps  are  possible  in  one,  two,  or  three  dimensions  and  allow  estimates  to 
be  made  on  size,  shape  and  motion  infornoation.  The  methods  vary  in  their 
use  of  data  with  certain  ones  depending  on  amplitude  or  an4>litude  and  phase, 
others  relying  on  polarization  or  frequency  diversity  and  some  requiring  a 
priori  motion  information.  However,  the  ability  to  perform  real  time 
imaging  so  necessary  in  the  present  applications  is  severely  restricted  or 
impossible  with  such  techniques  due  to  the  processing  time  requirements 
for  obtaining  the  various  single  dimension  of  resolution  and  the  need  for 
synthesizing  a  composite  multi-dimensional  map. 

Use  of  filled-in  large  apertures  constitute  the  basis  for  obtaining 
images  in  real-time.  Ilie  possibilify  for  obtaining  such  structured  three 
dimensional  images  using  radar-optics  analogs  directly  to  allow  immediate 
recognition  and  interpretation  establishes  the  necessity  of  such  an  approach 
in  the  realm  of  differentiation  and  identification  among  classes  of  target 
sources.  Of  course,  such  a  capability  can  also  provide  the  various  process¬ 
ing  alternatives  for  inverse  methods  already  mentioned  as  possible  with 
thinned  arrays. 
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Further,  the  use  of  microwave  methods  as  opposed  to  purely  optical 
methods  has  the  advantages  of  no  cloud  cover  and  greater  signal  to  noise 
sensitivity, 

FEASIBILITY  STUDY 

The  present  report  is  intended  to  provide  a  perspective  on  the  possible 
application  of  radar-optics  imaging  methods.  The  emphasis  wiU  be  on  those 
topics  which  have  been  treated  in  depth  to  date.  The  work  will  be  discussed 
in  relation  to  the  major  aspects  vlilch  justify  and  facilitate  the  application  of 
such  methods.  The  various  problem  areas  will  be  enumerated  and  relate  d, 
alternative  procedures  together  with  their  advantages  and  disadvantages  will 
be  discussed,  and  the  findings  in  particular  study  areas  will  be  presented  in 
sufficient  detail  to  establish  present  conclusions  and  future  efforts. 
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SECTION  n 
PRELIMINARIES 

BASIC  SYSTEM  CONSIDERATIONS 

To  be  successful  microwave  imaging  must  provide  a  display  which 
allows  visual  recognition  and  interpretation  of  the  scattering  source  in  at 
least  near  real-time.  It  is  useful  to  consider  target  sources  of  maximum 
characteristic  linear  dimensions  of  about  10  meters  and  to  provide,  for 
example,  a  possible  10:1  resolution  in  each  dimension.  The  object  to  be 
viewed  can  be  undergoing  translational  and  rotational  motion  either  alone 
or  in  the  presence  of  other  scatterers.  The  ranges  in  question  are  in  the 
order  of  a  few  hundred  kilometers. 

As  is  known,  for  such  ranges,  target  dimensions,  and  resolution 
requirements,  the  ground  array  overall  linear  dimensions  must  be  a  few 
tens  of  kilometers  and  the  array  stations  must  be  separated  by  no  more  than 
a  few  kilometers.  In  particular,  the  individual  array  stations  must 
be  spaced  sufficiently  clos^  so  as  to  sample  the  spatial  frequencies  in  vlie 
field  received  from  the  target  and  at  the  same  time  provide  a  sufficiently 
large  overall  array  size  to  obtain  the  desired  target  resolution. 

Of  course  array  thinning  can  occur  if  one  effects  resolution  by  a 
synthetic  aperture  approach  using,  for  example,  doppler  processing  at  the 
loss  of  real-time  operations.  It  is  the  desire  to  maintain  a  real  time  capa¬ 
bility  to  which  the  present  stu4y  of  a  large  real  aperatures  system  is  directed. 
Some  rearrangement  such  as  a  periodic  station  placement  may  be  useful  but 
for  the  present  report  only  uniform  spacing  will  be  considered. 
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The  scattering  at  mi  rowave  frequencies  expected  from  the  targets 
can  to  a  large  extent  be  adequately  represented  by  the  geometric  theory  of 
diffraction.  As  such,  the  target  is  basically  equivalent  to  a  collection  of 
scattering  centers.  Ideally,  and  assuming  a  sufficiently  large  dynamic  signal 
range,  the  number  of  centers  can  be  made  to  cover  the  body  by  means  of 
greater  diffusion  of  illumination  and  produce  an  image  approaching  that  ex¬ 
pected  at  optical  frequencies  at  least  to  the  extent  that  major  discontinuities 
on  the  body  will  appear.  To  what  extent  this  enhancement  of  image  quality 
can  be  attained  will  be  considered  as  part  of  the  continuing  investigation.  In 
any  case  the  different  bistatic  returns  obtained  by  each  array  station  from  a 
single  illuminating  source  can  cause,  for  example,  the  smearing  of  the  point 
scatterer  along  an  edge  and  so  provide  additional  visual  verificatton  of  an 
edge  geometry. 

The  application  of  holographic  methods  allows  the  possibility  of  direct 
three  dimensional  visual  displays  as  images  scaled  from  microwave  to  optical 
frequencies  without  the  need  for  explicit  processing  to  account  for  phase 
information.  This  requires  the  use  of  a  reference  field  in  addition  to  that 
from  the  target.  Such  p  reference  source  can  be  either  from  a  real  or  a 
synthetic  separate  source  or  possibly  one  derived  from  the  target  itself.  The 
former  approach  treated  in  the  present  report  establishes  strict  require¬ 
ments  on  the  knowledge  of  array  station  locations  to  small  fractions  of  the 
microwave  wavelength  as  well  as  the  need  for  a  fully  coherent  array.  How¬ 
ever,  it  appears  that  the  state  of  the  art  in  laser  ranging,  for  example,  can 
cope  with  the  station  location  requirement.  The  advantage  of  relaxing  know¬ 
ledge  of  station  location  accuracy  and  a  non-coherent  array  derive  from  use 
of  the  target  as  reference  source.  A  practical  scheme  to  accomplish  this  is 
still  lacking. 
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In  order  to  create  a  visual  image  from  the  collected  microwave  field, 
there  must  be,  as  noted,  a  conversion  from  microwave  to  optical  frequencies. 
Consideration  of  the  effects  of  using  sampled  data  to  create  a  continuous  view¬ 
able  image  is  to  be  given  large  emphasis  in  the  present  report.  Questions  on 
magnification,  resolution,  focus,  field  of  view,  and  image  displacement  are 
of  concern  here. 

In  order  to  practically  assess  the  various  i^stem  considerations,  an 
e;q}erimental  facility  capable  of  simulating  such  a  system  is  a  definite  require¬ 
ment  and  preliminary  investigations  on  such  a  capability  are  included  in  the 
study. 

MAJOR  STUDY  AREAS 

It  appears  that  a  practical  system  to  realize  the  desired  goal  of  ob¬ 
taining  improved  three  dimensional  images  in  real  time  can  be  approached. 

t 

This  possibility  is  based  partially  on  reasonable  estimates  on  array  size  and 
station  spacing  as  well  as  the  possibility  of  realizing  the  required  accuracy 
in  determining  station  location. 

To  appreciate  these  aspects  as  well  as  the  overall  investigation  more 
fully,  the  various  problem  areas  and  their  relationships  which  must  be  con¬ 
sidered  in  answering  questions  on  desirability  and  feasibility,  a  run  down 
of  the  major  study  areas  and  their  status  is  useful.  The  problem  areas  are 
grouped  generally  under  four  major  study  areas;  Justification,  feasibility, 
implementation,  and  interpretation. 

Justification 

A  large  portion  of  the  desirability  cf  such  an  imaging  method  hinges 
on  the  question  of  image  quality .  If  by  using  a  large  array  a  recognizable 
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visual  image  can  be  achieved  in  addition  to  having  a.  real-time  c^iabililyi 
then  such  an  approach  is  clearly  justified.  The  question  of  real-time  ciqia- 
bility  is  basically  established  by  the  large  array  arrangement.  Any  process¬ 
ing  times  associated  with  conversion  of  microwave  data  to  optical  frequencies 
should  not  appreciably  alter  this.  Actual  times  would  have  to  be  established 
after  the  details  of  processing  schemes  are  finalized. 

The  basic  question  of  image  quality  can  only  be  fully  considered 
combining  experimental  with  analytical  studies  which  must  be  in  part  approxi¬ 
mate.  The  need  for  establishing  an  experimental  capability  is  thus  para¬ 
mount.  In  fact,  the  lack  of  such  a  facility  to  date  has  prevented  a  fuller 
treatment  on  the  question  of  quality  in  the  present  report. 

Problem  areas  include: 

(a)  Scattering  range  work  (microwave  and/or  optical) 

range  design, 
experiment  design, 
simulation. 

(b)  Image  interpretation 

recognition  and  discrimination  of  shapes. 

The  choice  of  microwave  or  optical  simulation  e^qperiments  is  not 
clear-cut  and  will  be  discussed  in  the  section  of  the  present  report  titled 
E.xperimental  Work. 

Feasibility 

The  use  of  radar-optics  methods,  as  noted,  involves  the  use  of  a 
large  array  operating  at  microwave  frequencies.  The  process  of  forming 
an  image  after  conversion  to  optical  frequencies  can  be  accompP*^hed  as 
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either  a  one-step  or  two-step  process.  The  one-step  process  treats  the  array 
as  a  lens  while  the  two-step  process  is  the  holographic  method.  Considera¬ 
tion  of  both  methods  has  been  made  in  relation  to  questions  of  accuracy  of 
station  location  and  image  formation  and  resolution. 

The  present  report  considers  a  number  of  items  concerned  with  feasi¬ 
bility  such  as  the  various  aspects  of  one-or  two-step  array  imaging  methods 
and  a  comparison  of  these,  requirements  on  and  achievement  of  station  loca¬ 
tion  accuracy  atid  phase  coherency  across  the  array,  signal  sensitivity  in 
particular  with  respect  to  thermal  noise,  and  off- zenith  behavior.  Additional 
work  on  atmospheric  degradation  is  also  in  progress  but  will  be  detailed  in  a 
later  report.  The  question  of  a  practical  method  for  using  the  target  as  a 
reference  source  is  also  being  considered. 

Problem  areas  Include: 

(a)  Array  Imaging 

one-step  (lens), 
two-step  (holography). 

(b)  Station  Location  (array  geometry) 

requirements, 
feasible  methods. 

(c)  Coherent  Array 

requirements, 
feasible  methods, 
reference  signal. 


(d)  Atmospheric  Degradation 
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(e)  Sensitivities 

S/N  -  outputs, 
system  parameter  needs, 
coherent  integration. 


Implementation 

It  is  clear  that  a  number  of  items  must  bear  on  question'^  of  feasibility 
and  implementation.  Tlv^se  include,  for  example,  sensitivity  and  station 
location.  A  large  area  also  in  this  category  is  that  of  sanq>led  holography. 
The  present  report  denotes  considerable  attention  to  sampled  holognq)hy, 
giving  a  number  of  alternatives,  while  some  discussion  is  given  to  the  other 
items  in  the  list  below. 

Problem  areas  include: 

(a)  Sampled  holography 

method  implementation, 
image  effects. 

(b)  Sensitivity 

system  parameters. 

(c)  Other  Items 

aperiodic  array, 
reference  sources, 
illumination  diffusing, 
array  thinning. 

Interpretation 

The  question  of  output  usage  in  the  present  context  demands  that 
proper  interpretation  be  made  on  the  final  visual  image.  Work  presently 
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being  carried  on’*'  could  lead  to  an  approximate  analytic  basis  for  Judging  the 
effect  of  image  formation  accounting  for  polarization.  Possible  results  on 
these  matters  will  be  given  in  a  later  report.  Needless  to  say,  discrimination 
studies  given  under  Justification  also  bear  on  the  present  topic  of  interpretation. 
However,  this  and  the  other  items  on  the  list  below  will  be  dealt  with  in  future 
efforts. 

Problem  areas  include: 

(a)  Polarization 

image  enhancement. 

(b)  Optical  processing 

spatial  filtering, 
display  distortion. 

(c)  Computer  Processing 

filtering, 

reconstruction  quality, 
simulation. 

(d)  Techniques  for  improving  recognlzability 

smear  and  fill-in  effects. 

PLANS  FOR  FUTURE  WORK 

There  are  a  number  of  areas  which  are  to  be  treated  in  the  continuing 
work.  In  view  of  the  completed  investigations  and  present  status  our  first 
emphasis  will  be  to  undertake  the  investigation  of  image  quality  and  the  areas 
of  investigation  associated  with  setting  up  experimental  efforts  under  Justifi¬ 
cation  in  addition  to  work  in  the  other  categories.  This  will  include  questions 

♦By  R.  P.  Porter 
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on  work  in  station  locations:  non-planar  arrays,  generation  and  placement  of 
reference  sources,  array  thinning  and  aperiodic  spacing,  coherent  integration, 
chaff  effects,  diffuse  illumination,  polarization,  recognition,  and  discrimina¬ 
tion  processing  effects  on  real-time  cs^jability. 

In  addition  to  answering  questions  in  these  areas,  a  decision  on  one- 
step  or  two-step  processing  will  also  be  reached. 

REPORT  OUTUNE 

The  remainder  of  the  report  will  discuss  the  results  of  analytical 
studies  and  experimental  considerations  as  separate  items.  The  analytical 
work  will  center  around  three  main  topics:  array  imaging,  sampled  recon¬ 
struction,  and  sensitivity  and  S/N  effects.  These  are  discussed  in  summary 
fashion  in  the  body  of  the  report  followed  by  corresponding  detailed  analysis 
in  tne  Appendices.  Other  ^stem  considerations  are  given,  followed  by  a 
system  example  which  brings  together  all  of  the  analytical  work  and  provides 
a  set  of  system  parameters. 

The  next  discussion  in  the  body  of  the  report  concentrates  on  experi¬ 
mental  work  and  considers  requirements,  advantages,  choice  of  frequencies 
as  well  as  an  appropriate  optical  simulation  facility.  The  main  part  of  the 
report  ends  with  conclusions  based  on  completed  investigations  higddighting 
prospects,  deficiencies  and  future  directions  of  effort. 
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SECTION  ra 

SUMMARY  OF  ANALYTIC  STUDY 

RADAR-OPTICS  IMAGING* 

Introduction 

The  imaging  schemes  involve  a  basic  microwave  environment  con¬ 
sisting  of  an  array  surface,  S^,  situated  over  a  plane,  S^ ,  and  a  complex 
disturbance  located  in  a  finite  region  removed  from  S^  and  caused  by  the 
interaction  of  a  scattering  source  with  incident  quasi-monochromatic  radia¬ 
tion,  say,  of  wave  number  k  -  2ir/X  .  U  a  plane,  S^  ,  is  taken  parallel 

to  S  at  a  distance  i  above  S„  and  intersecting  the  region  of  complex 
H  H 

disturbance,  then  on  S  the  complex  distrubance  is  denoted  by  D(i)  .  Each 

o 

point  on  S^  is  given  by  vector  i  *■  and  on  S^^  by  vector 

X  =  (Xj^  ,  Xg)  .  ff  r'  is  a  vector  defining  a  position  on  the  surface  S^  , 
then  the  disturbance  on  S^  causes  a  field  on  S^  which  can  be  given  by 


where, 

1  is  a  unit  vector  in  the  z  direction  and 
z 

f ,  r' ,  S.  ,  and  S  etc.  are  shown  in  the  cross  section  that 
A  o 

follows: 


^Further  details  are  found  in  Appendix  1. 


12 


M68-2 


The  effect  of  the  disturbance  not  on  the  plane  can  be  treated  by 

considering  the  disturbance  on  a  sequence  of  such  planes  taken  at  different 

distances  from  S  given  by  f  +  .  in  general  ♦  (f '  ^))  *  ♦  (Jt)  and  for  S. 

XI  o  o  o  A 

planar,  ♦  (r')  =  ♦  (x) .  The  equation  for  ♦  reduces  to 
0  0  o 

i  .  r 

for  k »  1/r  and  -  *  1 

This  equation  is  the  basic  equation  which  is  the  starting  point  for 
either  a  one-step  (lens)  treatment  or  a  two-step  (holographic)  analysis.  A 
brief  indication  of  the  one-step  approach  to  illustrate  this  relationship  will 
now  be  given.  Since,  however,  the  bulk  of  the  stuc(y  in  this  reporting  period 
has  been  concerned  with  the  two-step  process,  the  major  portion  of  this 
summary  will  be  devoted  to  results  from  the  holographic  analysis. 
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Thus  a  simple  Fourier  pair  relationsh4>  exists  between  ^  and 

,2,, 

D(4)  If  the  phase  term  — jj—  ®  can  be  accurately  generated. 

V 

This  requirement  in  turn  demands  accurate  position  location  Icnowledge  in  x 
on  the  aperture  plane  to  a  fraction  of  a  wavelength.  When  the  Fourier  relation¬ 
ship  holds,  we  have  the  array  acting  as  a  simple  lens  equivalent.  Retrieval  of 
D(i)  visually  is  then  possible  by  simple  lens  Imaging  of  the  field  con¬ 

verted  to  optical  frequencies.  Such  a  process,  as  briefly  describee.,  is  the 
one-step  or  lens  imaging  approach  in  radar-optics. 


Two-Step  Process;  Hologram  Formation 


Returning  to  the  equation  for  'i'^(r*)  and  introducing  a  second  target 
as  a  point  source  reference  the  total  field  becomes 


» 


where  B  is  a  complex  constant  and  r^  is  the  distance  from  the  reference 

source  to  S . . 

A 

The  reference  is  located  at  (^_  ,  I,  _)  and  is  at  a  distance  I  +  from 

at  Jli  3R 

Then  without  imposing  the  Fraunhofer  condition  but  assuming  I  large 
in  comparison  to  the  array  and  object  extensions  in  the  phase  terms  and 
using  f  for  r  and  I  +  for  r-  in  the  denonoinators  of  i'  and  'i' 

on  M  OR 

respectively  the  formation  of  the  hologram  is  given  by 


IS 
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l.i 


(9) 


where,  using  r'  =  r*^)  and  denoting  I  as  I(x)  , 


X,cx)  Aj  ]  D^s)D  (f“)c  ^  Ax  Jtv 


'VS 


A^,L  LD  (S)« 


I.  (X)  sr  Aj 


(10) 


and 


A- 


JL 


,  .,  .;» cx«-  +  '*1^21  -  a.! 


a..bb!. 


(11) 
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This  step  of  constructing  the  hologram  together  with. a  second  step  of 
reconstructing  an  image  constitutes  the  two-step  or  holographic  imaging 
approach  in  radar-optics.  j 

Notice  that  with  the  Fraunhofer  condition  imposed  the  usual  Fcvier 
pair  relationship  exists  between  I  or  I  and  D(i)  if  compensation 
is  made  for  the  respective  phase  factors  A  or  A  =  A  ♦  . 

o  M 

Two-step  Process;  Reconstruction 

In  the  reconstruction  from  the  hologram,  it  is  the  satisfying  of  the 
focusing  condition  which  eliminates  the  effect  of  quadratic  phase  terms  and 
provides  effectively  a  Fourier  pair  relaticmship  and  a  resulting  image  of  the 
scattering  source.  In  dealing  with  a  synthesized  reference  source,  the  loca¬ 
tion  of  receivers  in  the  aperture  must  still  be  known  to  accuracies  of  a  small 
part  of  a  wavelength. 

The  reconstruction  will  take  place  at  optical  frequencies  so  that  a 
scaled  version  of  the  hologram  I^)  is  used  in  say  a  u  plane  as  F  (u)  with 

— A)  (12 


where 

m  is  the  scale  factor,  u^  a  bias  term,  and 
A  (x)  a  distortion  term. 

The  reconstruction  is  effected  by  parsing  a  spherical  coherent  wavefront  say 
of  wave  number  k'  through  the  scaled  hologram  F  (Q)  and  focusing  by  means 
of  a  lens  (with  focal  length  f).  This  leads  to  a  dl.'^tribution  4'(o)  in  an  image 
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plane,  ,  with  the  distance  from  the  lens  focal  point  to  the  scaled  holo 
gram  plane,  ,  and  the  distance  from  the  scaled  hologram  plane  to 
the  image  plane.  Using  range  approximations  with  respect  to  hologram  and 
image  extent  similar  to  those  in  the  microwave  situation,  the  image  field 
becomes 


-ih' E:^ 
e  »» 


JLSi 


(13) 


where 

<>'1; 

Q  -  e 

T  .  -iKAC* 

A  B  transmittance  ratio  from  hologram  (scaled) 

C  ■  intensity  of  reconstruction  light  source 
Cf 

-  >  intensity  of  light  impinging  on  hologram  (scaled) 


The  image  field  can  be  rewritten  for  the  separate  terms  in  I  as 


18 


e  r«  dJS, 


(14) 


M68-2 


With  i  =  2  for  example, 


Under  Uie  following  conditions. 


hlfl 


«  1 


Fraunhofer  condition  for  object 


^  ^  i-  2LpW  1  conditions  on  array  planarity. 

r  ^  i  L  / 


ik'Cf  *o 

X  (  (,  -f  5j^) 


optical  focusing  condition 


then  with  ^  =  0 


CS)s  Fx  ♦/m.TsaUw  \ 


where  F  '  is  a  complex  coefficient  in  which  a  only  appears  as  a  phase  factor 


so  that  (a)  is  magnified  by 


a 


mkz, 

T?r 


and  translated  by 


f  ^  « 


(on  the  object  plane)  image  of  the  original  object  distribu- 


3R 


tion  D(|)  .  The  image  is  real  for  z  >0  and  virtual  (viewable  image)  for 

V 

The  above  result  assumed  also  that  the  planes  and  (and  so 

^Hd ^  infinite  in  extent,  which  in  fact  violates  the  condition  e  i  . 

Taking  D(i)  =  0  outside  the  target  extent  causes  no  difficulty  but  since  the 
aperture  in  S„  (and  S„ .  )  is  of  finite  extent,  the  actual  value  for  I  should 
be  modified  by  a  truncation  function  in  the  x  plane,  say  Rect  (S)  which 
results  in  a  convolution  v  .th  a  sine  (^  function  in  ip.  (^  and  thus  a  "smear  out" 

h 

of  the  image.  Resolution  thus  improves  as  the  size  of  the  aperture  increases. 


Effects,  of  course,  from  ,  i  =  1,  3,  4  are  also  obtained.  For 


1  =  1, 


F>'(  £f(i)D(l-  - 


(17) 


which  is  a  convolution  smearing  of  D(i)  in  the  image  plane.  By  shifting 


f  ^  ^^3R  ^d 

(5)  off  axis  by  means  of  the  translational  term  -  ; -  (that 

2  Sr 

is  by  i-  or  u .  control)  then  a  non  overlap  between  and  is 

R  Q  ^2 

accomplished. 
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For  i  4  , 

<31  -  fv  s  (  ^  ) 


(18) 


a  delta  function  (ideally)  on  axis  in  the  image  plane.  Finally  for  i  =  3  , 


-f  ^  \ 


(19) 


where  with  focus,  (5)  is  not  (unless  =  0)  and 

Ip.  appear  on  opposite  sides  of  the  optical  axis  to  prevent  overlap. 
*3 


and 


Image  Characteristics 

The  various  conditions  for  the  selection  of  reference  placement, 
station  accuracy,  and  system  parameters  to  establish  an  image  in  terms  of 
the  microwave  receiving  aperture  array  and  the  optical  reconstruction  are 
illustrated  in  detail  in  the  section  on  a  Illustrative  System  Example. 


Using  the  basic  equations  and  conditions  noted  above,  various  para¬ 
meters  describing  the  imaging  process  can  be  derived.  Since  details  for 
some  of  these  are  given  in  Aj^ndix  1,  it  suffices  here  to  state  the  results. 
If  21.  is  the  overall  linear  dimension  extent  of  the  array  and  with  the 
Fraunhofer  condition  for  the  object  holding,  the  resolution  in  the  target 
plane  is  given  by 


Lj  =  Ai  , 
2  L 


(20) 


21 


and  in  the  z  (longitundinal  direction)  by 


J  j  ' 

L*  U 

which  is  much  poorer  than  L .. 

d 

Values  for  the  parallel  (to  the  object  plane)  and  the  perpendicular  (in 
the  z  direction)  magnification  factors,  which  are  ratios  of  distances  in  the 
image  plane  over  respective  distances  in  the  object  plane,  can  be  given  as 

,i  -5  >>«■ 

"  '“FT 


and 


s  />r»„ 

Kr-  ~l~ 


(22) 


SAMPLED  PROCESSING  AND  HOLOGRAPHY* 

To  obtain  microwave  information  by  an  array  aperture  for  eventual 
conversion  to  an  optical  image,  a  sampling  of  the  received  field  is  required. 
In  order  to  accoimt  properly  for  the  spatial  frequencies  corresponding  to  a 
target  source  of  approximately  10  meters  in  extent,  a  minimum  spacing 

* 

Further  details  are  found  in  Appendix  2 
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distance,  (2b),  is  needed,  and  should  be  S4)proxlinately  a  few  kilometers. 

The  spatial  extent  cf  each  receiving  aperture,  ustially  less  titan  a  few  meters, 
for  example,  must  be  small  compared  to  the  minimum  spatial  period  caused 
by  any  increase  in  effective  ao«mce  extent  resulting  from  a  combination  of 
target  and  reference  source.  Finally,  the  number  of  receivers  is  dett-r- 
mincd  by  the  overall  extent  of  the  array.  For  resolution  cells  in  the  order 
of  1  meter,  this  requires  array  dimensions  of  a  few  tens  of  kilometers  or 
approximately  100  sampling  stations  over  the  array  surface. 

Once  the  field  is  accurately  measured  in  amplitude  and  phase,  various 
possibilities  can  be  used  to  apply  this  sampled  data  to  form  a  continuous 
visual  image.  The  first  and  simplest  method  so  far  considered  and  discussed 
in  detail,  converts  in  a  straightforward  manner  from  a  sampled  microwave 
hologram  to  a  sampled  hologram  scaled  to  optical  frequencies.  Over  each 
sampling  ceil  in  two  dimensions,  the  hologram  amplitude-only  function,  I, 
is  fixed  at  its  center  value.  Reconstruction  takes  place  directly  from  the 
sampled  hologram. 

Two  alternatives  have  also  been  analyzed,  and  constitute  the  second  and 
third  methods.  The  second  method  is  more  complicated  than  the  first  since 
it  converts  the  sampled  hologram  mentioned  above  into  a  continuous  hologram 
as  an  intermediate  step,  prior  to  reconstruction.  The  third  approach  (sxig- 
gested  by  R.  Roig)  first  converts  the  sampled  field  in  amplitude  and  phase 
into  a  continuous  field  prior  to  any  processing  for  a  visual  image.  This 
method  has  the  major  disadvantage  that  control  of  both  the  amplitude  and 
phase  of  the  very  small  reconstructing  light  sources  must  be  maintained 
over  the  scaled  aperture.  More  particularly,  extremely  severe  requirements 
are  placed  on  the  location  of  these  li^t  sources  as  scaled  down  from  the 
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microwave  location  requirements.  However  sampling  amplitude  and  phase 
of  the  field  rather  than  the  intensity  used  i;i  holography  has  the  advantage  of 
requiring  the  minimum  number  of  samples  without  introducing  error. 


To  begin,  the  first  method  converts  the  hologram  I(x)  sampled  at 
positions  (x  ,  x  )  in  x  on  S  into  the  scaled  version  (scaling  factor  m), 

1]  II 

with  sampling  cells  centered  at  (u,.  ,  u.  )  in  u  on  S„.  where,  with  p, 

1]  2n  nd 

and  A  -  0,  Xjj  .  m  Ujj  and  -  m  . 

The  sampling  takes  place  uniformly  at  distance  2b  in  each  dimension 
of  X  ,  and  I  is  taken  constant  over  each  sampling  cell  so  that  for  all  J,  n 
within  the  cell  defined  by: 

where  b'  may  be  less  than  or  equal  to  b, 

X  (7)  *  I  (X.i, 

with 


Thus,  in  the  S„  plane,  the  sampled  version  of  I  results  in 
H  2 


and,  in  the  S  plane, 
Hd 


x,'  (0=  'W 

••  It  ^ 
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The  reconstructed  field  then  bficomes,  in  the  Sj  plane, 


»iy 


a  (26) 


where 


b^ 

Together  with  l“i  -  “ijl  "=  .  ^“2  "  ^  5  ’  to  permit  focusing, 

the  following  conditions  must  be  satisfied. 


where 


( 

^  * 


J. 


Taking 


4 


and 


U,M,  ‘t 


(28) 


then  the  reconstructed  image  becomes 
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T»,j*/»w*-  e  •** 


X2  I(«.«,1,-H..V 


If  the  conditions 


(30) 

are  assumed,  and  D(|)  is  taken  to  be  zero  for  1{  ^  1  >  ,  i  ■  1,2, 

where  2Lj^  is  the  maximum  linear  extent  of  the  target  in  each  dimension  of  the 
of  the  plane,  then  in  the  region  (cj  ^  * 

D(f)=  (3» 


0  (sU. 

Sli  l,  ^  ^ ' 


where  D  is  ^e  Fourier  transform  of  D. 

On  inserting  the  appropriate  focusing  condition,  and  taking  b  to 
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then 


*2^)=  YM  -  -fijl 


(33) 


a).-k*B*ACiX 


■ailA 


(34) 


which,  with  the  various  assumptions  and  restrictions,  such  as  on  b,  dem¬ 
onstrates  the  possibility  of  forming  a  magnified  and  translated  image  of 
D(|).  The  modulation  intensity  due  to  the  use  of  a  sampling  cell  hologram 
is  clearly  indicated  in  the  sine  functions.  To  reduce  degradation,  the  image 
must  be  contained  between  zeroes  of  these  fimctions  in  the  a  plane.  Of 
course,  the  always-present  smearing  effect  of  the  finite  array  size  will 
cause  additional  degradation,  as  was  pointed  out  in  the  previous  subsection 
on  radar-optics  imaging. 

The  transition  to  the  above  expression  for  (a)  from  the  proceeding 

o2 

expression  for  that  quantity  is  presented  in  detail  in  Appendix  2;  however, 
it  is  appropriate  to  point  out  here  that,  althou^  it  is  not  apparent  from  the 
above,  the  fields  in  the  image  plane  resulting  from  the  various  terms  of  the 
hologram  are  periodic,  so  that  fields  from  the  different  terms  interfere  with 
each  other.  This  interference  can  be  reduced  by  focusing. 
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The  question  of  degradation  of  the  image  represented  the  above 
expression  wUl  bo  dealt  with  at  greater  length  in  Appendix  2  . 

In  addition  to  the  sampling  requirement  stated  above,  the  conditions 
I  2  2|  ,2  2  I 

on  1  Uj^  -  u^j  I  and  [  previously  given  imply  restrictions  on 

|^.„|  for  a  given  value  of  b.  This  restriction  may  be  stated  approximately 
««  2Ij  i 

(assuming  I «  I).  as  . 

It  is  assumed  in  the  foregoing  that  the  scattered  field  from  the  target 
(and  from  the  reference  source,  if  an  actual  reference  is  used)  is  sampled 
at  a  point.  This  requires  that  the  dimensions  of  the  receiving  antenna  be 
small  compared  to  the  reciprocal  of  the  highest  spatial  frequency  appearing 
in  the  measured  field  (i.  e. ,  small  compared  to  ,  where  is  the 

maximum  linear  dimension  of  the  complex  of  target  and  reference,  if  there 
is  an  actual  reference. )  This  is  not  a  severe  restriction. 

Additional  considerations  on  the  significance  of  the  various  approxi¬ 
mations  is  required,  for  example,  in  any  comparison  of  the  various 
sampling  methods.  Some  discussion  on  the  requirements  of  power  levels, 
coherent  integration,  noise  levels,  target  and  reference  source  constraints, 
and  image  size,  is  given  in  the  sections  on  Sensitivity,  Further  System 
Considerations,  and  Illustrative  System  Example. 

Alternative  Approaches 

In  order  to  avoid  the  sine  modulation  and  the  periodicity  in  the  image 
plane  caused  by  the  finite  sampling  cell  size,  the  second  reconstruction 
method  first  converts  the  sampled  hologram  to  a  continuous  hologram.  This 
is  done  by  first  taking  the  Fourier  transform  with  an  "f  to  f"  lens  for  each 
sampled  hologram  cell.  The  new  field  is  then  given  by 
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"'▼‘^  w  . .,  •i->«i2ddsHMiw 


^  Irtrv  — 
(36) 


If  ihe  focal  plane  of  these  processing  lenses  is  such  as  to  let  zeroes 

of  sine  function  occur  at  w„  =>  u. ,  ,  then 

J’n  j*n’ 


(36) 


and  the  sampling  theorem 


J/**- 


X(»> 


(37) 


that  is,  a  field  having  magnitude  prc^rtional  to  the  original  continuous 
hologram  (scaled  down)  is  produced; 


.  iik!f 


(38) 


Use  of  this  field  permits  reproduction  of  a  continuous  hologram  from 
which  an  image  of  the  original  distribution,  D(| )  may  be  obtained. 

The  third  method  seeks  to  r^roduce  the  (scaled  down)  original 

scattered  field,  ,  directly.  This  leads  to  using  tho  set  of  values 

^  (X  )  as  weii^ts  for  the  interpolation  with  the  effect  of  a  multi-hole 
J I  n 

pierced  lens  to  generate  sine  interpolation  functions,  as  were  used  in  the 
second  method,  to  reproduce  now  ^(x) .  The  complication  of  implementing 
and  locating  a  complex  light  modulation  at  the  holes  with  values  ^  ^  ) 

j 

is  most  severe. 
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T1ii6  approach  Involves  smearing  effects  fxom  finite  cell  size  and 
depth  of  focus  difficulty.  An  implementation  to  facilitate  obtaining  a 
viewable  image  has  been  suggested.  Details  on  such  matters  will  not 
be  presented  in  tiie  present  report. 

SENSITIVITY  AND  THERMAL  NOISE* 

Of  particular  concern  in  studying  signal  degradation  is  the  effect  of 
receiver  thermal  noise  on  the  image  reconstructed  from  sampled  data.  This 
can  be  measured  in  terms  of  a  signal-to-noise  ratio  (SNR)  between  image  and 
noise  intensity.  It  appears  likely  that  this  ratio  should  be  largely  unaffected 
by  the  particular  method  used  for  image  reconstruction.  The  effective  SNR 
depends  in  a  complex  manner  on  the  nature  of  the  scattering  source  for 
given  array  and  receiver  noise  characteristics.  This  ratio  of  image  inten¬ 
sity  to  background  receiver  thermal  noise  intensity  is  strongly  influenced  by 
the  nonlinearities  of  the  holographic  process. 

In  order  to  expedite  obtaining  an  appreciation  of  the  importance  of 
SNR  and  its  influence  on  system  parameters,  a  number  of  approximations 
are  made.  The  method  used  is  generally  applicable  and  gives  results  in 
terms  of  such  system  parameters  as  transmitter  power,  transxnitter  and 
receiver  antenns.  gains,  array  dimensions,  and  coherent  integration  times. 

In  an  attempt  to  investigate  the  effects  over  a  range  of  targets, 
approximate  expressions  for  the  SNR  are  given  for  an  unresolved  point 
target  and  a  resolved  plate  target.  A  function  of  the  reference  signal  level, 
thin  ratio  asymptotically  ajqproaches  a  maximum  as  the  reference  level 
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increases.  For  the  unresolved  target,  the  ratio  depends  <m  die  soattering 
cross  section,  while  for  the  resolved  plate,  the  asymptotic  value  is  inde¬ 
pendent  of  plate  size. 

Ihe  Noise  Model 

The  analysis  begins  with  a  number  of  assumptions.  All  antennae  in 

the  array  receive  the  same  polarization  under  matched  load  conditions. 

The  effective  received  voltage  comes  from  N  coherently  added  pulses. 

Ihe  resultant  noise  is  essentially  unvarying  in  time  intervals  less  than 

s  1/A  f,  where  the  receiver  bandwidth,  Af ,  is  small  compared  to 

center  frequency  f  ,  and  r  is  single  pulse  duration. 

P 

This  leads  to  received  voltage  given  by 

where  J  refers  to  the  particular  pulse,  V,.  (»  for  all  j  ,  assuming 
appropriate  phase  shifts  to  permit  coherent  integration)  is  the  received 
pulse  complex  voltage  without  thermal  noise  and  the  received  complex 
noise  voltage  for  the  Jth  pulse  is  n^  =»  ^  ^nj^yWith  jn^j  and 

and  <l>  essentially  constant  during  the  Jth  pulse. 


The  ensemble  average  on  lv^|^  gives 

(V»f  =  N*lV-rl‘+ 


where 


/rij 


(independent  of  J) 
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The  voltage,  NV^  ,  would  be  produced  by  a  single  pulse  incident  electric 
field  of  magnitude 


(41) 


where  R  is  the  radiation  resistance  of  the  antenna,  G  the  gain  of  the 

A 

receiving  antenna,  and  tj  the  impedance  of  free  space. 

^  Also,  the  incident  electric  field  corresponding' in  the  same  sense, 
N 

to  the  noise  voltage,  ^  n^  ,  1ms  mean  square  mag^tude  of 

J-1 


.  tirhUlT^i , 

(S» 


(42) 


where  k  is  the  Boltzmann  constant  and  T  is  the  total  noise  temperature. 
N 

Thus  the  voltage,  ^  n^  ,  present  after  coherent  integration  of  N  pulses, 
J-1 

can  be  jrepresented  for  narrowband  thermal  noise  by 


where  the  amplitude,  >  is  Rayleigh  distributed^ with  |e^|  given 

above,  and  the  phase,  0  ,  is  uniformly  distributed  over  (o ,  2ir ),  with 
(Ej^I  and  ^  independent  at  different  stations  and  independent  of  each  other 
at  a  given  station  for  thermal  background  noise. 
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The  effective  total  field  at  the  array  plane  la 

E/7)  +  E^<S)  +  EnCt),  (43) 

(object^  (/n**^*) 

and  the  power  flux  la  |E^l^/2n,  with  j  a  measure  of  the 

hologram  intensity ,  I .  By  proper  adjustment  to  separate  the  terms  of  the 

hologram  reconstructed  image  can  be 

obtained  from  the  term  EE*.  Now,  however,  the  following  degradation 

D  K 

terms  caused  by  the  noises  appear  in  the  hologram. 

'  E„l*-  +  E,  E-*-*  E^e!  ♦  E,  Ef  +  £h  eJ  (44) 


Noise  Effect  on  Reconstructed  Image 

In  the  analysis  of  the  reconstruction  process  firom  sampled  data, 
the  expression  for  i/i  (a)  contained  a  sum,  S,  which  has  a  mean  squared 

o2 

magnitude  of 


(45) 


With  assiimptlons  of  hoinogenity  and  independence,  an  evaluation  results  in 
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S  S*  =  m(i  Ewl'j*  +  2  lE.l*'  2  P 


where 


\n  “  ®pj,n  ®RJ.n 


M  total  number  of  stations 

1  z  *  ^■•***‘  !*■ 

Ji<K 

Considering  P  as  a  random  walk  near  the  edges  gives  P  on  the  order  of 
M  .  Thus, 


Assuming  additional  reduction  in 
over  ]  ,  n  gives 


due  to  phase  variation 


J  .  n 


S  5*  =  M 1  e«l‘  C  X 16/^'  ;  +  2  (e J*  +  E Ja  ) j 
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where 

With  no  thermal  noise  of  course,  S  S*  ■  0  . 

Application  to  Two  Target  Sources 

Two  specific  targets  are  now  considered:  (1)  Point  target  unresolved, 
and  (2)  flat  plate  resolved  target. 

Point  Target  Unresolved 


It  is  assumed  that  the  target  is  located  at  L  with  D(|) 

s 

A'  6(i  -  Ig)  ,  A'  =  constant  and  the  array  surface  planar  with  x.  »  0. 
The  target  field  (f ' )  *  E^(f ' )  ,  and  target  cross  section  is  cr .  A' 
is  evaluated  from  energy  flow  considerations  as 


where 

P^  -  transmitted  power  during  each  pulse, 
>  transmitter  gain, 

I  ^  =  range  from  target  to  transmitter. 

After  N  coherent  returns. 


(51) 


(52) 
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The  reference  field  at  the  array  ^  (r*)  =  E  (?')  is  given  by 

K  R 


E,«)=  B  LIk£_l 

*  *‘*-*^>  A+5»  »(i+T,0 


With  B  constant. 


I  1 2  2 

With  I E  I  and  |  E  j  assumed  constant  over  the  array 
»e  Rj ,n 

(i.e.  over  j  ,  n),  giving 

the  mean  squared  magnitude  of  the  contribution  to  the  reconstructed 
field,  (a)  ,  from  the  contribution  in  S  due  to  the  noise  term  in  the 
hologram  1 .  is  given  by 

1E.«  1"=  Sr.KS?  ,5s 


where  the  sine  terms  in  ^  are  neglected  since  in  any  SNR  calculation  they 

IS2 

would  not  appear. 

To  obtain  a  SNR  a  comparison  is  made,  after  reconstruction,  between 
the  squared  field  magnitude  of  the  center  of  the  diffraction  pattern  of  the 
holographic  aperture  (which  is  actually  obtained  in  reconstruction  rather 
than  a  point  image)  with  the  mean  sqiuired  magnitude  of  the  noise  field. 
Using,  for  a  large  number  of  stations,  the  continuous  hologram  in  place  of 
the  sampled  hologram  to  reconstruct,  becomes 
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E  ^  ^  V  (5e) 

where 

1F.|*=  (A-)\^iBi^(i>;)M^ei" 

i=.^  *>V|  —  k<*»5’M  \  ni7\ 

^  jTtn"  J 


('  _  km.'%u.\ 

The  squared  reconstructed  field  at  the  point  where  the  image  of  the 
point  target  should  appear,  after  N  coherent  pulses,  is 

|E„,  f-  L»lBr(hrA*C‘-(-*N'-Pr67<rn 

a  Kiif’  U*  S,*)’- 

The  SNR  is  then 


(59) 
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Flat  Plate  Resolved  Target 


It  is  assumed  that  D(0  *  B'  Rect.  (4/a)  and  the  perfectly-conducting 
plate  has  sides  of  length  2a  >>  X  so  that  the  broadside  cross  section  is 
64ira^ 


From  power  flow  considerations 


and  after  N  coherent  returns  with  a  Fraimhofer  conditon  on  the  target, 


The  array  average  ,  ,  can  be  ai^roximated  a  continuous 

average  over  the  area  of  the  array  to  yield  Ejj^^  ~  N^P^G,pha^/2Til 
over  the  array,  then  neglecting  the  effects  of  sine  terms,  as  witii  the  noise 
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As  iBf  -  «>  , 


(depends  on  a) 


K. 


'h*i'FK 


(independent  of  a)  (65) 


where 

^  Qr<^iL  N 

Note  that  the  magnitude  of  the  reconstructed  field  at  the  center  of  the 
diffraction  spot  for  a  point  target  equals  the  magnitude  of  the  reconstructed 
field  in  the  interior  of  the  image  for  a  plate,  if  the  point  target  has  cross 


section 


or-  Th'-jj 
-f  u- 


2  2 

This  also  gives  R.  =  R_  .  Of  course,  these  values  can  in  general 
UMAX  i-<MAX 

differ  widely  as  the  numerical  example  in  Aiq}endix  3  clearly  indicates.  For 
example,  if  N  =  100 


-2  2 

with  a  =  10  m  for  the  point  target  and  the  resiilt  being  independent  of 
the  plate  cross-section,  provided  that  the  plate  is  large  enou^  to  be 
resolved.  Such  a  spread  in  SNR  suggests  the  need  for  some  type  of  ampli¬ 
tude  compression.  Additional  computation  also  indicates  conqutibility  between 
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total  transmit  time  for  N  =  lOO  pulses  and  the  time  to  restrict  relative 
movement  between  major  target  scattering  points  to  be  less  than  l/io  of 
the  operating  wavelength. 
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SECTION  IV 

FURTHER  SYSTEM  CONSIDERATIONS 

STATION  LOCATION  REQUIREMENT  FOR  LARGE  APERTURE  RADAR 
ARRAYS 

The  analysis  of  radar-optics  imaging,  as  summarized  in  Section  m 
and  detailed  for  holographic  images  in  Appendix  1,  indicates  that  in  obtaining 
the  image^unless  an  actual  reference  source  is  used,  one  of  the  most  crucial 
requirements  imposed  on  the  system  for  one-  or  two-step  imaging  is  the  need 
of  knowing  the  location  of  each  station  to  within  a  small  fraction  of  the  opera¬ 
ting  wavelength. 

Once  the  location  of  all  stations  is  known  to  such  accuracy,  it  is  then 
possible  to  synthesize  a  reference  near  the  target  by  tracking  the  target  and 
calculating  the  appropriate  complex  reference  field  that  each  station  would 
have  observed  if  the  reference  point  source  were  actually  so  located. 

Althou^  the  above  requirement  seems  quite  stringent,  it  is  neverthe¬ 
less  possible  with  the  present  technological  means  to  obtain  the  location  of 
each  station  within  the  desired  accuracy.  This  report  will  briefly  review 
long- baseline  measuring  techniques^  ,  the  station  location  accuracy  re¬ 
quirement  and  then,  to  establish  feasibility,  will  indicate  a  systematic 
method  of  measuring  the  accurate  location  of  a  large  number  of  stations. 

Accurate  Distance  Measuring  Techniques 

There  are  two  types  of  long-baseline  measuring  techniques:  those 
using  the  "time  of  arrival"  ai^roach  and  those  using  the  "phase  difference" 
approach;  furthermore,  both  of  the  above  types  can  be  subdivided  into  two 
broad  classes:  optical  and  microwave  measuring  schemes. 
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The  "time  of  arrival"  technique  at  optical  or  microwave  frequencies 
measures  distances  by  sending  a  short  pulse  signal  to  a  planar  reflector,  lo¬ 
cated  at  one  end  of  the  path  to  be  measured,  and  recording  the  transit  time 
back  to  the  transmitting  or  a  separate  receiver,  located  at  the  other  end  of 
the  path. 

The  geometrical  path  D  between  the  above  two  points  is 

_  r  ^ 

(68) 

where  t  is  the  two-way  transit  time  and  V  is  the  average  group  velocity  of 

8 

the  transmitted  short  pulse  over  the  path.  The  range  accuracy  depends  on 
the  accuracy  of  measuring  the  transit  time  and  the  average  group  velocity  V  . 

D 

Although  the  accuracy  in  transit  time  gets  better  as  the  pulse  length 
gets  smaller,  a  too  short  pulse  will  be  broadened  by  a  dispersive  medium 
such  as  the  atmosphere  and  hence  there  is  a  limit  on  how  much  can  be  gained 
in  accuracy  by  going  to  a  smaller  pulse  length. 

The  accuracy  of  optical  range  measurements  by  means  of  pulsed  lasers 
is  better  than  that  of  microwave  measurements  because  of  reduced  atmos¬ 
pheric  effects  due  to  a  lesser  water  vapor  contribution. 

Before  discussing  atmospheric  effects  and  other  systematic  errors  in 
general,  let  us  indicate  the  other  approach  which  is  used  to  measure  precisely 
long  distances. 

The  "phase  difference"  or  "Interferometry"  technique  was  first  investi- 

[21 

gated  by  Bergstrand  in  the  early  1950*  s'-  ,  and  it  laid  the  foundation  of  modern 

optical  distance  measurements.  The  approach  can  briefly  be  summarized  as 
follows:  a  light  beam  is  amplitude-modulated  at  a  frequency  f  and  sent  out 
to  an  optical  reflector  at  the  end  of  the  path  to  be  measured. 

The  returned  light's  phase  is  compared  to  that  of  the  transmitted  light. 
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If  the  measured  phase  difference  is  Ad ,  the  two-way  path  lengtti  to 

the  end  reflector  is  2L  =  (n  +  X  and  it  is  ambiguous  because  of  the 
unknown  quantity  "n"  . 

The  ambiguity  can  be  removed  by  performing  the  same  measurement 
at  various  frequencies. 

A  similar  scheme  can  be  implemented  for  the  microwave  system. 

Very  precise  geometrical  path  distances  between  two  points  can  only  be 
obtained  if  all  forms  of  errors  are  properly  accounted  for.  The  factors 
affecting  the  accuracy  of  measurements  for  both  optical  and  microwa\'e  sys¬ 
tems  are: 

(a)  Stability  of  transmitting  source, 

(b)  atmospheric  effects  due  to  the  changing  index  of  refraction, 

(c)  accuracy  of  the  phase  measuring  system,  and 

(d)  uncertainty  from  lunar  and  solar  tides  in  the  solid  earth. 

For  lung  distance  measurements,  the  atmospheric  effect  is  the  biggest 
problem. 

Taking,  as  an  example,  the  "phase  difference"  approach,  the  path  L 
that  is  obtained  is  the  optical  path  which  is  the  sum  of  two  terms:  the  geo¬ 
metrical  path  D  plus  the  additional  path  difference  Q  due  to  the  atmosphere. 

L  =  D  Q 


L 


i)  as 


(69) 
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where  n'  the  group  index  of  refraction  defined  as 


V,(s) 


is  a  varying  function  along  the  path.  To  obtain  the  desired  geometric  path, 
it  is  necessary  to  calculate  the  optical  correction 


which  can  only  be  done  if  the  group  index  of  refraction  is  known  at  many 
points  along  the  path.  However,  since  the  atmospheric  contribution  to  the 
index  of  refraction  is  about  300  parts  per  million,  the  air  density  must  be 
known  to  one  part  in  300  in  order  to  permit  an  overall  distance  accuracy  of 
a  part  per  million.  Thus  the  atmospheric  pressure  and  temperature  for 
many  points  along  the  path  must  be  known  to  one  part  in  300,  a  rather  diffi¬ 
cult  task.  To  avoid  making  temperature  and  pressure  measurements  at 
many  points  along  the  measured  path  in  calculating  Q  (the  optical  correc- 

[3  4] 

tlon)  a  scheme  has  been  recently  proposed  for  the  optical  case^  '  ^  . 

It  uses  measurements  obtained  at  two  largely  different  wavelengths 
(one  in  the  red  and  the  other  in  the  blue  part  of  the  visible  spectrum);  the 
approach  takes  the  difference 


s  Q,  -  Q, 

of  the  two  optical  path  corrections  (equal  to  the  optical  path  difference  AL) 
and  finds  that  the  ratio  of  AQ  to  is  a  constant  ivhich  is  essentially 
independent  of  the  atmospheric  density  along  the  path  and  is  known  from 
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fc  g1 

laboratory  measurements''  *  ^ .  Thus  an  accurate  m^isured  value  of  AQ  can 
be  used  to  find  the  desired  path  correction  to  approximately  the  same 
fractional  accuracy. 

As  an  example,  for  a  path  of  15- km  AQ  is  about  40  cm,  and  this 

7 

would  have  to  be  measured  to  0.4  cm  to  give  3  part-in-10  accuracy 
in  L 


AQ  _  ,  -> 

-  =  ■  .  — 2.  s  - - “7  arlitu 

L  AQ 

A  similar  scheme  could  be  devised  for  the  microwave  system.  How¬ 
ever,  because  of  the  necessity  of  using  two  greatly  different  wavelengths  and 
therefore  two  transmitting  and  receiving  systems,  its  practicability  is 
greatly  reduced. 

In  addition  to  the  above  drawbacks  the  microwave  distance  measuring 
systems  seem  to  be  less  accurate  than  the  optical  systems  because  of  their 
greater  sensitivity  to  atmospheric  water  vapor. 

The  overall  accuracy  of  present  optical  long-baseline  distance  measur- 

0 

ing  systems,  correcting  for  the  atmospheric  effects,  is  one  part-ln-10  . 

If  the  scheme  using  two  greatly  different  wavelengths  is  used,  then  it 

7 

is  hoped  that  an  overall  accuracy  of  a  few  parts  in  10  can  be  achieved. 

Station  Location  Accuracy  Requirements 

When  a  large  radar  array  is  used  to  do  optical  imaging,  there  are  two 
modes  in  which  such  array  can  be  operated:  either  as  a  lens,  or  as  a  holo¬ 
graphic  recording  medium.  In  both  cases,  however,  the  requirement  on 
sution  location  can  be  obtained  from  the  following  simple  analysis. 
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Given  an  object  field  distribution  D(i)  a  distance  I  away  from  the 
recording  non-planar  array  and  assuming  for  the  time  being  a  continuous  array 
(see  diagram  below),  the  complex  field  produced  at  point  P  is  (See  Section  m 
and  Appendix  1) 


I 
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Assuming 


liili" 


«1  ftnd 


JU.  I 
Jf 


-  -.*• 


Using  the  above  ai^roximation  in  the  e3qx>nent  and  the  approximation  r  »  I 
in  the  denominator,  the  integral  becomes 


'®  iirf 


If  a  proper  phase  compensation  is  introduced  across  the  array  such  that  the 
new  complex  field  amplitude  over  the  array  is 


,  -‘ir  liTKivfjfrt) 

t  =  J  ^3><-  ^  ^  Jv  , 


and  if  the  Fraunhofer  condition  is  applied  Lil  44  1  .  together 

2.  f 


with  the  array  planarity  condition  : 


x'l  (?•  l) 
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then 


(X)  ^  J  «-  *  *% 


(79) 


and  a  Fourier  transform  condition  exists  between  the  complex  field  ip^  and 
the  object  field  distribution  D(|)  . 

In  effect,  however,  because  of  station  location  inaccuracy,  the  phase 
term  over  the  array  will  deviate  from  the  desired  one  and  will  add  a  random 
phase  contribution  to  the  useful  complex  field  . 

Such  random  phase  distribution  will  degrade  the  quality  of  the  Fourier 
transform  of  and  thus  will  decrease  the  resolution  of  the  image. 

A  criterion  must  then  be  imposed  on  the  accuracy  of  measuring  the 
location  of  all  stations  which  clearly  depend  on  how  much  phase  variations 
one  can  tolerate  over  the  array. 

Looking  at  the  phase  term 

if  (8») 

e  =  € 


the  phase  error  is  clearly  related  to  the  errors  in  the  measurements  of  f  ,  X 

and  X 
3 
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The  maximum  phase  error  Is  obtained  when  x  is  maxinuim: 

^max  ~  ^  (where  2L  is  the  array  size)  and  when  all  errors  are  arranged 
so  as  to  give  such  a  maximum, 

S  k[A»  +  f  ~ 

since  *V  ^  ^ 

A*  :s  W.[Af  +  .1.  t  AXjl 

MM  L  ^ 


However,  in  determining  L,  one  needs  the  vertical  position  of  all 

AX3 

M  stations  of  the  array  and  therefore  Af  »  ~5l“  •  For  M  “  100, 
the  M  term  can  be  neglected  and 

AA  =  K  ^  (83) 

From  the  above  relation  one  sees  that  distance  errors  in  the  x  -  direction 
are  ten  times  more  critical  than  those  in  the  x  -direction  and  contribute 
mostly  to  the  phase  errors. 


49 


M68-2 


AssunUng  that  one  can  measure  Ax  and  Alxl  with  the  same  accuracy, 

then 


k 

In  the  following  analysis  we  assume  that  a  maximum  phase  error  of 
20^  will  not  affect  the  quality  of  the  reconstructed  image.  In  the  future,  a 
more  detailed  investigation,  considering  the  random  phase  errors  over  the 
array,  will  be  carried  out  to  determine  exactly  what  maximum  phase  error 
one  can  tolerate  in  order  to  be  able  to  successfully  reconstruct  the  image. 
For  the  moment,  however,  taking  A<b  at  20°,  the  maximum  tolerated 
error  in  distance  is 

AX*  ^  -ii.  -  -L  >s  *5  ..5  for 

Actually,  the  distance  error  in  the  planar  direction  needs  to  be  about  ten 
times  less  critical 


2r  10 

The  fractional  accuracy  in  the  planar  direction,  if  L  w  15  km,  is  then 

Afft  .  5  _  _ 

L  i.sxu* 

an  accuracy  which,  using  the  optical  distance  measuring  system  outlined 
above,  can  presently  be  achieved. 
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A  Systematic  At 
Radar  Array 


iroach  for  Locating  Large  Number  of  Stations  In  a 


This  approach  will  briefly  outline  a  method,  using  the  accurate  dis¬ 
tance  measuring  Instruments  discussed  previously,  to  measure  the  location  of 
all  stations  of  the  radar  array  to  within  the  desired  accuracy. 


For  a  raaor  array  of  100  stations  or  more,  the  task  of  precisely 
locating  in  three  dimensions  all  stations  is  a  long  and  involved  process. 


It  is  desirable  to  keep  the  number  of  optical  distance  measuring  Instru¬ 
ments  to  a  minimum  and  to  set  14)  the  e^qperiment  so  that  the  errors  in  the 
planar  direction  are  about  as  equal  as  those  in  the  vertical  direction. 

Bouncing  laser  beams  from  station  to  station  will  give  a  good  accuracy 
in  the  planar  direction  but  not  in  the  non-planar  or  vertical  direction. 


On  the  other  hand,  bouncing  laser  beams  ofi  satellites  will  decrease 
the  accuracy  in  both  dimensions  because  of  the  long  paths  measured. 

A  compromise  which  is  proposed  here  is  to  bounce  laser  beams  be¬ 
tween  stations  themselves  as  well  as  between  stations  and  nonrplanarly 
located  targets  (such  as  taU  towers,  balloons  or  helicopters);  this  ^^roach 
will  maintain  the  accuracy  desired  and  in  addition  will  equalize  the  errors 
in  the  different  dimensions. 

Furthermore,  to  keep  the  number  of  optical  instruments  to  a  minimum, 
it  is  useful  to  locate  a  minimum  number  of  stations  at  a  time.  First,  the 
whole  array  can  be  divided  in  several  subgroups,  each  made  up  of  such  a 
minimum  number.  After  each  station  of  a  subgroiq)  is  located  relative  to  a 
"common"  station  of  that  subgroup  for  all  subgroups,  the  next  step  would  be 
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to  locate  in  position  and  orientation  each  subgroup  relative  to  a  "common" 
subgroup. 

Let  us  then  determine  the  minimum  number  of  stations  which  are 
sufficient  to  permit  determination  of  their  locations  relative  to  one  another. 

Each  station  has  an  optical  distance  measuring  instrument  bouncing  a 
laser  beam  off  reflectors  fixed  on  the  other  stations  and  off  reflectors  fixed 
on  a  non-planarly  located  target  (balloon,  tower,  helicopter,  etc. ). 

The  independent  unknowns  are  the  3-dimensional  locations  of  all  N 
stations  relative  to  one  of  them,  i.  e. ,  a  "common"  station,  and  the  3-dimenr 
sional  location  of  the  non-planarly  located  reflector  on  the  target. 


(85) 


The  number  of  different  measurements  that  can  be  obtained  bouncing  laser 
beams  off  reflectors  fixed  on  the  other  stations  and  the  non-planarly  located 
target  is 

IP  ^  ^ 


n(n-i)  ^  ^  •  tiOil!) 

Z  2 


(86) 


To  be  able  to  solve  for  the  unknowns,  the  following  inequality  must  be  imposed 
between  the  number  of  unknowns  and  the  number  of  measurements. 
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or 

iKUlO  V  2  fN/ 

or 

M  ^  s 

Thus,  using  one  measurement  per  station  and  a  non-planarly  located 

target,  at  least  5  stations  are  required  to  completely  determine  their  position. 

For  a  large  array  of  M  stations,  using  the  above  approach  and  using  five 

M 

laser  distance  measuring  instruments  -r  times,  all  stations  on  each  of  the 
M  ® 

~  subgroiq)s  will  be  located  relative  to  one  "common"  station  of  each  sub- 
5 

group. 

The  next  step  is  to  locate  in  position  and  orientation  all  subgroups 
relative  to  a  "common"  subgrotgi.  This  can  be  done  first  by  locating  in 
position  any  orientation  sets  of  subgroups  relative  to  a  "common"  subgroup 
of  each  set,  and  then  all  of  the  sets  relative  to  a  "conunon"  set. 

Since  we  are  locating  in  orientation,  as  well  as  in  position,  three 
points  of  each  subgrotg)  must  be  located  each  time  ^’elative  to  a  station  of 
the  "common"  subgroiq). 

Conclusions 

This  report  has  considered,  ;n  some  detail,  station  location  accuracy 
requirements  for  optical  imaging,  the  various  distance  measuring  techniques 
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one  could  use  to  achieve  such  requirements,  and  finally  has  outlined  a  system¬ 
atic  approach  of  getting  accurate  three-dimensional  locations  of  a  large  num¬ 
ber  of  stations  using  a  minimum  number  of  q;)tical  distance  measuring  instru¬ 
ments.  There  is  further  work  to  be  done  in  investigating  how  the  random 
phase  errors,  introduced  over  the  array,  decrease  the  resolution  of  the 
reconstructed  image  and  how  the  mayhinnTn  station  location  accuracy, 
affecting  the  resolution,  limits  the  maximum  resolution  of  the  radar  array  for 
coherent  optical  imaging,  and  puts  a  limit  on  the  use  of  the  array,  beyond 
which  no  gain  in  resolution  is  attainable. 

OTHER  SYSTEM  CONSIDERATIONS 

There  are  many  other  system  considerations,  which  although  quite 
important,  have  not  as  yet  been  investigated  in  detail. 

Nevertheless,  it  is  felt  that  some  of  these  considerations,  which  will 
contribute  additional  areas  of  research,  should  be  discussed  here  to  indicate 
possibilities  of  sinq;)lifying  or  improving  the  imaging  ability  of  the  discussed 
large  radar  array. 

The  four  main  considerations  which  will  be  described  briefly  here  are: 

(a)  Coherence  Requirements  and  Possibility  of  Achievement 

(b)  Reference  Source  Used:  Synthetic  or  Actual 

(c)  Diffuse  Illumination 

(d)  Off-Zenith  Behavior  of  the  System 
Coherence  Requirement  and  Possibility  of  Achievement 

As  shown  in  Section  m,  if  a  large  radar  array  is  used  as  a  lens  to  do 
optical  Imaging,  the  basic  requirements  imposed  on  the  system  are  a  precise 
knowledge  of  station  location  and  a  ftilly  phase-coherent  array. 
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The  theoretical  analysis  developed  in  the  subsection  dealing  with  station 
location  accuracy  has  shown  that  any  inaccuracy  in  the  determination  of 
station  location  will  add  a  random  jdiase  variation  over  the  array  which  will 
degrade  the  reconstructed  image.  Similarly,  any  inaccuracy  in  obtaining  a 
fully  coherent  array  will  also  add  a  random  phase  over  the  array  and  thus 
will  also  affect  the  resolution  of  the  reconstructed  image.  The  total  random 
phase  variation  over  the  array  considered  as  a  siun  of  the 

two  above  contributions: 


(88) 


Assxuning  that  the  maximum  allowable  rms  total  phase  error  at  each  station 
is  20°,  and  if  the  phase  error  due  to  station  location  inaccuracy  is  10°,  all 
stations  must  be  phase-coherent  to  within  10°.  Again,  for  a  large  number  of 
stations,  this  is  a  quite  stringent  requirement  which  nevertheless  can  be 
achieved  in  several  different  ways. 


One  approach,  which  is  customarily  used  to  provide  a  phase  reference 
for  multistatic  radar  systems,  is  the  use  of  phase  stable  microwave  links. 
For  distances  greater  than  a  few  kilometers,  fluctuations  in  the  transmitting 
mcdiiun  necessitate  multipath  procediu'es  to  correct  for  such  fluctuations. 


A  better  approach  of  obtaining  a  coherent  phase  reference  is  by 
utilizing  the  frequency  stability  of  atomic  clocks.  These  are  accurate  to 
1  part  in  10^^.  If,  for  each  station  of  the  radar  array,  an  atomic  clock  is 
provided  and  all  clocks  are  calibrated  in  phase  by  observation  of  an  astro¬ 
nomic  source,  a  coherent  phase  reference  would  be  established  over  the 
array. 

An  interferometer  using  atomic  clocks  to  provide  coherent  local 
oscillator  signals  has  been  operated  for  astronomical  purposes  witb  a  base- 
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f  7  1 

line  of  461, 000  X  at  610  MHz  ^  ^  yielding  resolution  only  sli^tly  inferior 

f  ft  1 

to  that  attainable  with  optical  systems  ''  ^  . 

One  proposal  considers  the  construction  of  a  3000-mile  baseline 
interferometer,  using  atomic  clocks  to  provide  a  phase  reference.  At  the 
suggested  operating  frequency  of  1000  MHz,  phase  recalibration  would  be 

required  about  once  a  day.  With  the  same  clock  stability,  operation  at 
3000  MHz  (S-band),  recalibration  would  be  required  about  three  times  as 
often.  The  clock  accuracy  is  more  than  adequate  for  our  purpose  and 
when  using  svich  an  approach  the  total  phase  error  will  mostly  be  due  to 
station  location  Inaccuracy. 

A  third  approach  woxild  be  a  combination  of  the  two  previously  discussed 
methods.  It  would  consist  of  noaintalning  phase  coherence  of  several  subgroups 
of  stations  by  microwave  links  and  phase  coherence  between  each  subgroup  of 
the  array  atomic  clocks.  All  the  different  approaches  will  be  investigated 
more  thoroufdily  in  future  work. 

Reference  Source  Used;  Synthetic  or  Actual 

When  the  radar  array  is  used  as  holographic  recording  system,  the 
requirements  of  station  location  accuracy  and  coherence  accuracy  depend 
on  whether  an  actual  or  a  synthesized  reference  point  source  is  used. 

If  an  actual  reference  point  source  is  used  near  the  target  and  is 
moving  with  the  target,  each  station  directly  measures  the  amplitude  and 
phase  difference  needed  for  obtaining  the  hologram  and  thus  a  coherent  phase 
reference  is  not  required  over  the  array.  Each  station  must  still  be  capable 
of  measuring  the  phase  and  anmlitude  of  the  complex  field.  However,  the 
system  needed  to  preserve  coherence  throufidxout  the  array  is  eliminated. 
Furthermore,  since  each  station  now  measures  (mly  the  relative  phase 
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between  the  object  and  the  reference  point  and  not  the  absolute  phase  of  the 
object  scattered  field,  and  since  such  relative  phase  depends  to  a  lesser 
extent  on  station  location  if  the  reference  is  close  to  the  target,  the  require¬ 
ment  of  station  location  is  also  relaxed.  To  investigate  how  much  this 
requirement  is  relaxed,  let  us  consider  the  complex  object  field  e}q)ression 
over  the  radar  array  obtained,  when  the  Fraunhofer  condition  and  the  array 
planarity  condition  are  obeyed. 


i 

XTTjj 


I 


t  < 
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(89) 


If  a  reference  point  is  present  at  the  origin  of  | -plane,  the  complex  field  of 
that  reference  la  ^ 

-B  .  -B  -i _ ^ 

r 


%  - 


■«  ( 

The  virtual  image  term  (see  Section  III  and  Appendix  I)  can  be  written 


as 


fpW  -  -  -LL-  i  1>(J)  C  * 

iT  J 


(91) 


and  it  is  directly  related  to  the  Fourier  transform  of  the  object  field  dis¬ 
tribution  D(i)  .  Required  station  location  accuracy  depends  now  only  on 
the  variation  of  the  above  term  and  not  on  an  additional  phase  term  in  front 
of  the  integral  (Equation  89).  The  maximam  variation  of  the  above  equa¬ 
tion  or  the  smallest  fringe  over  the  array  occurs  for  the  largest  target  dimen¬ 
sion  ^  max. 
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For  the  one-dimensional  case  and  assuming  D(|)  =  const  =  K  over 
the  target,  the  above  integral  becomes 


-  /OO' 


— 

2. 

The  lobe  width  Ax^  of  which  is  obtained  by  setting 

!lJ^  AV.  a  T 

2  ^ 


7  3 

I  «  300  km  -  3  X  10  cm  {  =  10  cm  X  =  10  cm 

max 

.  ikla^xie  »  J  T  to 

^  - —  s  i  XIo  cw«  =  3 


— 

Choosing  one  one  hundredth  of  such  width  as  the  criterion  for  station  location 
tolerance,  one  still  sees  that  the  requiremmit  becomes 


AX  « 


t  V  w 

clearly  a  much  more  relaxed  requirement  than  if  an  actwil  reference  was  not 
near  the  target.  The  actual  reference  from  a  target  could  be  a  prominent 
scattering  center  or  a  specular  point  on  the  target.  However,  in  most  cases, 
an  actual  reference  is  hard  to  find  because  a  specular  point  mig^t  not  be 
present  for  that  view  and  a  prondinent  scattering  center  might  not  be  visible 
at  all  stations  of  die  array. 
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If  a  synthesized  reference  is  used  to  obtain  a  hologram  (or  for  that 
matter,  if  the  array  is  used  as  a  lens)  then  it  must  be  possible  to  calculate 
the  phase  variation  over  the  array  due  to  a  point  source  located  near  the 
target.  The  phase  variation,  given  by  Equation^(90)  is 

^  +  ]  m 

and  it  changes  the  field  values  function  at  the  array  to  a  new  field 

value  function  (x)  ,  which  is  the  Fourier  transform  of  the  object  source 
distribution  D(i)  .  (See  Equation  78). 

,  f  ,  -“li  - 

t.  (>■)  =  J  (1.)  ‘ 

To  be  able  to  calculate  such  a  phase  variation  over  die  array,  the  stringent 
station  location  and  coherence  requirements  of  the  one-step  processing  must 
be  imposed  on  the  system. 

Diffuse  Illumination 

When  optical-type  imaging  is  done  for  general  targets  by  means  of 
large  aperture  radar  arrays  operating  at  microwave  frequencies,  the  image 
obtained  resembles  that  which  would  be  obtained  from  very  smoodi  objects; 
more  specifically,  only  a  few  highli^ts  (such  as  specular  points,  edge 
scattering  centers)  show  iq)  on  die  object,  while  the  remaining  features 
remain  hidden. 

The  reason  for  such  peculiar  behavior  of  microwave  imaging  versus 
optical  Imaging  for  the  same  target  is  that  the  microwave  wavelength  used 
is  much  larger  than  the  optical  one,  and  thus  most  information  is  obtained 
from  regions,  aside  from  specular  points,  which  have  discontinuities  or 
rou^mess  larger  than  the  microwave  wavelength  used. 


For  one  transmitter  and  a  particular  receiver,  the  number,  position 
and  type  of  scattering  centers  on  a  general  target,  which  contribute  mostly 
to  the  electromagnetic  scattering,  and  thus  those  points  which  will  image, 

[  9  1 

are  given  by  Keller's  geometrical  theory  of  diffraction 

Using  an  array  of  receivers,  but  still  only  one  transmitter,  each 
receiver,  seeing  the  target  from  a  slightly  different  direction,  will  measure 
the  scattering  due  mostly  to  sli^tly  different  points  on  the  target.  Hence, 

the  image  obtained  will  be  a  more  filled  one  than  for  monostatic  systems, 
since  each  scattering  center  will  be  smeared  along  the  edges  or  disconti¬ 
nuities  on  the  target,  and  thus  will  add  to  the  recognizability  of  the  target. 

In  order  for  the  microwave  image  to  resemble  the  optical  image  of 
an  object,  it  is  necessary  to  simulate  the  diffuse  scattering  at  optical 
wavelengths. 

One  approach  toward  that  objective  consists  of  illuminating  the  target 
fkvm  different  directions,  so  that  although  for  each  illuminator  the  target 
does  not  scatter  diffusively,  the  combined  scattering  for  all  illuminators 
a{^)ears  diffuse  to  the  extent  that  returns  from  the  target  come  from  the 
large  fraction  of  the  extent  of  each  region  of  discontinuity,  and  not  only  from 
few  scattering  centers. 

If  an  image  is  then  processed  usir^  the  array  data,  the  image  will  look 
like  an  outlined  sketch  of  the  target,  resembling  the  target  more  closely 
than  would  a  three-dimensional  structure  of  scattering  centers. 

There  are  several  ways  of  implementing  such  "diffuse"  target 
illumination.  Using  about  a  dozen  transmitters  placed  at  different  locations 
(some  of  which  will  be  outside  the  array)  one  could  illuminate  the  target 
simultaneously  or  in  sequence. 
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The  simultaneous  illumination  requires  that  the  different  transmittere 
produce  a  combined  incident  wavefront  at  the  target,  which  is  diffuse  but 
coherent  at  the  same  time  in  order  to  be  able  to  do  coheroit  itnegiiig 
using  the  radar  array. 

This  requirement  imposes,  in  turn,  stringent  station  location  accuracy 
as  well  as  coherency  requirements,  so  that  appropriate  phase  delays  may  be 
used  between  transmitters  to  permit  the  construction  of  the  incident  coherent 
wavefront.  Furthermore,  since  the  targets  dealt  with  by  the  array  are 
moving,  an  extremely  precise  tracking  of  the  target  is  required  in  order 
to  vary  the  phase  delays  between  transmitters  in  order  to  always  be  able  to 
construct  a  coherent  incident  wavefront.  One  of  the  many  difficulties  of 
generating  such  diffuse  illumination,  using  the  simultaneous  approach,  is  the 
additional  phase  change  introduced  by  atmospheric  inhomogeneities  along  the 
different  paths  for  each  illuminating  beam,  degrading  the  combined  coherent 
wavefront.  A  quantitive  analysis  of  the  degradation  of  the  coherent  wave- 
front  generated  by  several  transmitters  due  to  atmospheric  inhomogeneities, 
together  with  other  considerations  involved  in  producing  a  coherent  diffuse 
illumination,  will  be  conducted  in  the  future. 

The  sequential  illumination  approach,  not  having  to  generate  a  combined 
coherent  incident  wavefront,  does  not  suffer  from  the  station  location  and 
coherency  requirements,  but  on  the  other  hand  does  involve  the  fact  that 
since  each  transmitter  illuminates  in  sequence  the  moving  target,  it  will 
be  viewed  at  different  times. 

The  question  of  maximum  coherent  integration  times,  pulse  duration, 
dealt  with  in  Appendix  3,  must  enter  in  the  discussion  of  the  sequential 
approach  to  simulate  diffuse  illumination. 
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For  target  velocity  V  «  7000  m/«eo,  the  pulse  duration  time,  t  , 
necessary  for  the  target  not  to  have  moved  more  than  1/10  X  (X  »  10  cm) 
in  that  time  interval,  was  found  to  be 

T  »  '0 

The  xnaximum  coherent  integration  time  T  ,  to  allow  for  as  hl{^ 

max 

a  S/N  ratio  as  possible  witibout  blurring  the  hologram,  was  found  to  be 


based  on  estimates  of  maximum  target  size  and  rotation  rate.  The  slgnal- 
to-nolse  analysis  indicated  the  desirability  of  coherently  integrating  at  least 
100  pulses. 

Using  the  sequential  iqpproaoh  to  diffuse  illumination  without  having 
to  maintain  phase  coherence  beMreen  stations,  the  most  convenient  way  to 
proceed  mlg^  be  to  transmit  from  each  transmitting  station,  for  example, 
one  hundred  1  p-second  pulses,  coherently  integrate  them  at  the  receivers, 
process  them  to  obtain  the  hologram  and  reconstruct  an  imag".  then  super- 
inrpose  10  of  these  separate  images  to  obtain  tiie  "combined*'  final  image  which 
will  be  a  more  detailed  sketch  of  the  target  than  would  each  of  the  individual 
images. 

A  blurring,  due  to  the  motion  of  the  target,  mi|^t  be  present.  To 
reduce  the  blurring,  one  mic^t  obtain  each  separate  image  wldi  less  coherent 
integration  time. 

All  of  the  above  approaches  to  sim\ilate  diffuse  illumination  will  be 
investigated  in  greater  detail  in  the  future. 


The  reeulte  obtained  in  the  teohnloal  aiuUysis  developed  sc  far  consider 
the  target  to  be  located  at  tiie  zenith  of  the  radar  array. 


In  most  oases,  however,  tiie  target  will  not  be  at  the  zenith  of  the 
array's  center  and  a  brief  analysis  is  needed  to  indicate  die  differences 
that  do  arise. 

Consider  the  same  diagram  as  included  at  the  beginning  of  this  section, 
but  with  a  target  placed  off-zenith  and  assume  for  simplicity  a  planar  arr^r. 
(See  diagram  below). 
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The  complex  field  ^^(x)  at  ai^  point  P  of  the  array  due  to  a  complex 
object  source  distribution  D(|)  located  in  the  vicinity  of  is  given  by 


the  usual  e}q)ression 


t.  -  - 


Y  -  - 

j 


where 


r 

K  r.  II 


(102) 


-  a 


Assuming  that 


i  .  .1 


where  i  'I*  «  I, 


the  above  expression  for  r  >  ftfter  a  change  of  variable,  ^  ^  > 


becomes 


J  !£L  - 
•  z  f, 
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Using  such  an  expression  for  r  in  the  e;q)onent  of  the  integral,  but  only 
the  approximation  r  in  the  denominator,  the  complex  field  ampli¬ 

tude  at  point  P  becomes 


*'  *’■  JUVsOt  ■' 

^  JTT  4,  -y 
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Assuming  the  Fraunhofer  condition, 

and  introducing  the  proper  phase  modulation  over  the  array,  one  can  obtain 

the  new  complex  field  ^  '  (x) : 

o 


(105) 


which  is 


(106) 


The  Fourier  Traneform  of  over  an  Infinite  planar  array  wIU  then  bet 
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which  is  the  complex  source  distribution,  D(i^  +  a)  ,  apart  from  a  linear 
phase  term.  Since  one  is  interested  only  in  the  magnitude  square  of  the  image 
and  not  in  the  complex  quantity,  the  phase  terms  will  not  effect  the  quality  of 
the  image. 

If  the  planar  array  is  of  finite  extent,  2L  ,  in  one  of  the  two  planar 
directors,  according  to  the  above  equation,  the  complex  source  distribution 
D($)  is  then  convolved  with  the  function 

tine  (1«« 


(108) 


The  transverse  resolution,  Lj  ,  in  the  cr-plane  for  a  point  source  in  the 

e 

$ -plane  is  then  obtained  by  setting 


(109) 
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The  above  e}q)ression  Indicates  the  the  off-zenlth  transverse  resolu¬ 
tion  is  poorer  than  the  zenith  transverse  resolution  by  a  factor  cos  0  , 
where  6  is  the  angular  orientation  of  the  viewed  target  relative  to  the 
zenith  of  the  array.  This  result  is  consistent  with  vae  optical  resolution 
formula,  which  equates  the  resolution  to  \/a  ,  where  a  is  the  angle  sub¬ 
tended  by  the  optical  viewing  system.  Similarly,  tlie  longitudinal  resolution 
i  ^  ,  and  the  parallel  and  perpendicular  magnifications  given  in  Section  m 
and  in  Appendix  1,  will  change  and  will  be  poorer  when  the  target  is  in  an 

off-zenith  position.  A  complete  investigation  of  such  off-zenith  behavior  will 
be  carried  out  in  the  near  future. 

An  additional  result  will  be  affected  by  the  off- zenith  behavior;  for  a 
target  located  at  an  angle  6  relative  to  the  normal  to  the  array,  the  sampling 
interval,  established  for  the  d  -  0  direction,  will  be  able  to  measure  higher 
spatial  frequencies  because  the  sampling  interval  will  appear  smaller  by  a 
fraction  cos  d  in  the  d-dlrection  and  thus  targets  of  larger  spatial  extent 
can  be  dealt  with.  In  summary,  the  off-zenith  behavior  of  the  system  will 
permit  imaging  larger  regions  of  space  with  poorer  resolution  and  magnifica¬ 
tion  at  low  elevation  angles,  and  imaging  smaller  regions  with  better  resolu¬ 
tion  and  magnification  at  hi^er  elevation  angles. 
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SECTION  V 

an  illustrative  system  example 

In  order  to  aUow  greater  appreciation  of  the  various  constraints  and 
requirements  associated  with  microwave  imaging,  calculations  for  appropriate 
parameters  are  given  that  are  applicable  to  a  two-step  process  with  sampled 
data  and  optical  reconstruction. 

To  begin  we  choose  the  following  operating  data: 

microwave  wavelength,  x  =  o.lm  (k  =  ^  =  62.8m" ^); 
optical  wavelength,  x  •  =  0.5  x  10"®m  (k*  =  — ,  =  1.26  x  lO^m"^; 
frequency  ratio,  k'/k  =  X/X*  =  2  x  10®; 
f  =  3  X  10®m; 

Ljj  -  5m  (half  linear  target  extent); 

=  Im  (linear  resolution  size). 
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MICROWAVE  SET-UP 


Array  Size  (2L): 

From  L,  = 
d 

2L  = 


30  km. 


Microwave  Station  i^aclng  (2b); 


ai»  i  -  3  w. 
2L^ 


Number  of  Stations  In  Square  Array  (M); 


ail) 


(112) 


ai3) 


M  »  ( 


2  L 
2  b 


12) 


(114) 


Obliquity  Condition  (Microwave); 

,  (L+L„t 
(A-xtY-  'jf~ 


(116) 


Fraunhofer  Condition  (Microwave); 


(116) 
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Toleraxice  on  Array  Planarity: 

Using  l~  I  +  4^,  and  assuming  |  <  1km. 


(  T»a  .  zLtN  ^  -  1 

r  '<T-  "  T^)  '  ^ 

requires 


\X^\  «  2  k'm. 


(117) 


(118) 


M6e>2 


/ 


so  it  is  required  that 


>r  ^ 


(122) 


In  which  |£._|/l  «  1  has  been  used.  Herce,  if  ^ =  “500m  , 

3xi  «ixi 


^123) 


or 


(124) 

/•n. 


7 

where  m  is  the  scale  factor,  here  set  equal  to  0.6  x  10  ,  because,  as  will 
be  seen,  such  a  value  leads  to  convenient  optical  parameters.  A  value  of  b' 
of  20m  would  be  reasonable.  A  condition  similar  to  the  above  pertains 
to  lug  -  I  ,  resulting  in  the  same  requirement  on  b*. 


Optical  Focusing  Condition 

The  focusing  condition  for 

S2  “ 


(t2S) 
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in  which,  it  will  be  recalled,  and  are,  respectively,  the  distances 
from  the  focal  point  of  the  spherical  wave  used  in  the  reconstruction  to  the 
hologram  and  the  distance  from  the  hologram  to  the  image.  A  negative 
value  of  z  corresponds  to  a  virtual  image.  Using  =  -500m, 

2^5  (126) 

If  Zj^  *  0.9m, 


(127) 


Image  Magniflcation  from  Reconstruction  and  Viewing 


^  0 

Overall  Magnification 


1/2  anrie  on  eve  (aided  or  not)  from  imafce 
1/2  an^e  on  eye  ^rom  original  object 


Let  /3  »  1/12  radian  «  5° 

6 


»  X“  K/O* 

^•/jt 


(128) 


But  m^  may  be  written  as  (Mag^)  (Mag^) 


and 


Mag-  the  anguls^r  magnifications  due  to  the  reconstruction  process  and  due  to 
viewing  the  reconstructed  virtual  image  with  a  telescope,  respectively.  A 
convenient  set  of  optical  parameters  results  from  a  choice  of  Mag^  -  30 
and  Magg  =  6  x  10^/30  »  166.4  . 
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Reconstruction: 


From  m— 
k' 


30,  m  =  0.6  X  10 


7 


(scaling  factor),  so  that 


A*'* 


(129) 


The  expression  mp 


for  Mag^  comes  from 


Ik* 

since  for  1  unit  change  in  a ,  image  plane,  there  results  ; -  units 

kmz 

change  in  |  ,  object  plane. 


Thus, 


i*  *  .  X 

‘-hx 


(131) 


OM-i 

CSCAl»i) 


73 


M68>2 


where  1 1  i  the  parallel  mr  ’ficatlon,  is  discussed  in  Appendix  1.  The 
angular  magnification,  Mag^,  is  then  obtained  as  follows.  Let 


a  is. 

X 


I'M 


(132) 


Then 


i-h  lta\  Kf 


(133) 


Viewing; 

Telescope  viewing:  objective  lens:  JL  4 


lU  U  +. 


eyepiece  lens:  *  —  4«  -J 


(135) 


With  telescope,  i  =  -  «  ,  so  that  i  =  f  , 
'  ®I  ®o  ® 


(136) 


A  Mag  of  166.4  may  be  achieved  with  Zg  ~  "9^1  =  6.8mm. 
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Non-Overlap  of  Images; 

For  non-overlap  of  1.  and  I  we  consider  the  following; 

X  £t 


Separating  the  I  and  I  images  will  automatically  separate  the  I.  image 

X  2  « 

from  the  and  images.  The  effective  length  of  the  image  from  is 
4L^m||  ,  since  it  is  the  autocorrelation  of  D(i) . 

If  we  let  Y  -  be  the  angular  displacement  of  the  center  of  the  image 

O 

due  to  from  the  zero  axis,  then  from  the  sketch  it  is  required  that 

Y,  >  ^  »  '5*  •  (137) 

From  the  form  of  ,  the  origin  of  coordinates  in  the  object  plane  is 
IS2 

mapped  into  the  point  a  in  the  image  plane  (considering  only  the  image 

02 

resulting  from  I^)  where 

77^  -  *  Pw) 
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ms» 


•WM-r 


so  that 


I  ^*>^1  „  ^AK  |Tr.\ 

I*.-!  ~  i  K' 


(139) 


assuming,  as  before  that  Ugj^l  «  f  .  Accounting  for  the  viewing  tele¬ 
scope,  we  want  an  angular  shift  such  that 


I  ®Cox\  Ij 

i*u.r 


t  4x 

Jl  4e. 


/wi -r  2  3 


(140) 


or 


I  A 


(141) 


Image  Placement  with  Respect  to  Zeroes  of  Sinx/x  Intensity  Modulation 


It  is  shown  in  Appendix  2,  Method  1,  that  the  reconstructed  image  is 
modulated  in  intensity  by  a  two-dimensional  sinx/x  bmction.  The  zeros  of 
this  function  are  separated  by 


or,  in  angle  by. 


(  AoC  I  C 

K'V 


l£il_  =  Cjl  i  S  wtfc  ^  fc 

l?b  *>  IT  “jr  T 
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in  which  the  sampling  condition  2b  =  lir/L^k  has  been  used.  It  will  be 
recalled  that  mk/k*  is  just  the  angular  magnification  resulting  ih:om  the 
reconstruction  process,  and  2L^/l  =  2|3^  ,  the  angle  subtended  the 
target  as  seen  by  an  observer's  unaided  eye.  Thus 


\^>C\ 


a44) 


where  2^  „  is  the  angle  subtended  by  the  virtual  image  of  the  target  as 
seen  by  an  observer  in  the  plane  of  the  scaled  down  hologram.  Thus,  if 
b/b*  2:  image  fits  between  zeros  of  the  intensity  modulation.  In 

this  example  b'  was  chosen  previously  to  be  *  20m,  and  b  was  chosen 

3 

as  1.5  X  10  m.  Hence,  b/b'  »  75,  and  the  variation  of  intensity  modula¬ 
tion  across  the  image  is  small. 
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EXPERIMENTAL  WORK 

REQUIREMENTS  AND  ADVANTAGES 

Experimental  measurements  are  required  in  order  to  investigate  and 
establish  actual  results  at  microwave  frequencies  that  are  independent  of 
any  analyses  which  rely  on  approximations.  This  need  is  most  critical  to 
supplement  the  analysis  that  provides  data  to  process  an  image,  judge  its 
quality,  and  evaluate  recognition  capability.  In  addition  of  course,  study 
and  verification  of  other  areas  associated  with  one-  or  two-step  processing 
can  be  carried  out.  Illustrative  of  the  wide  range  of  application  for  experiment 
are  type  of  illumination,  type  of  reference  and  placement  effects,  station 
location  accuracy,  array  thinning,  array  non-planarity,  sampling,  image 
reconstruction,  and,  generally,  evaluation  of  error  effects. 

Areas  that  could  be  investigated  concerning  image  quality  include  the 
kind  of  image  possible,  image  realism,  filled-in  view  of  target  discontin¬ 
uities,  monostatic  and  multistatic  illumination  with  simple  or  multiple 
transmission,  and  the  benefit  of  polarization  diversity.  In  addition,  methods 
of  carrying  out  coherent  processing  for  image  improvement  within  time 
constraints  compatible  with  target  movement  and  real-time  outyuts  are  of 
concern.  Chaff  effects,  image  distortion  and  discrimination,  and  the  use  of 
spatial  filtering,  for  example,  are  also  important  in  determining  inoage 
quality. 

Once  accurate  and  appropriate  microwave  data  is  available,  the  task 
of  synthesizing  such  data  and  evaluating  the  results  of  the  many  experiments 
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possible  can  be  handled  ly  computer  processing.  However,  even  though 
obtaining  experimental  data  is  fundamental,  questions  on  range  versatility, 
feasibility,  and  applicability  are  important. 

The  degree  to  which  optical  simulation  is  used  will  depend  on  (1)  whether 
or  not  range  data  is  taken  initially  at  optical  frequencies,  and  (2)  whether  one- 
or  two-step  image  processing  is  being  considered.  The  use  of  direct  optical 
processing  for  constructing  an  image  with  the  one-step  approach  is  not  as 
obvious  as  in  two-step  processing,  and  will  be  dealt  with  in  future  study, 
particularly  since,  as  noted  in  the  conclusions,  the  one-step  processing 
might  be  favored  over  the  two-step  processing,  when  all  factors  have  been 
considered. 

CHOICE  OF  RANGE  FREQUENCY 

The  desired  goal  is  to  perform  range  measurements  that  simulate 
the  actual  target-array  microwave  environment  as  closely  as  possible. 

Ihus,  it  is  desirable  that  the  array  be  in  the  target  far  field,  that  is, 
that  the  Fraunhofer  condition  hold. 

If  one  chooses  an  operating  range  frequency  in  the  microwave  spectrum 
under  such  conditions,  the  target-to-array  distance  and  array  size  become 
large  and  impractical.  Alternatively,  for  reasonably  sized  arrays  and  dis¬ 
tances,  operation  in  the  microwave  region  actually  requires  Fresnel  scatter¬ 
ing  rather  than  Fraunhofer  scattering. 

On  the  other  hand,  operating  in  the  optical  frequency  range,  for  which 
a  wealth  of  technical  and  experimental  data  exist,  permits  a  proper  Fraunhofer 
simulation  together  with  reasonable  and,  in  fact,  very  convenient  target-to- 
array  distance  and  array  size.  A  possible  difficulty  occurs  in  the  scaled- 
target  size  which  can  be  in  the  order  of  hundredths  of  a  millimeter.  Even 
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here,  a  large  variety  of  target  specimens  are  available  to  represent  actual 
targets  ranging  from  smooth  ones  to  those  with  many  discontinuities. 

Considering  the  factors  of  overall  convenience  and  proper  simulation, 
it  appears  that  experimental  data  on  the  target  can  be  taken  conveniently 
at  optical  frequencies. 

OPTICAL  SIMULATION 

Once  the  data  is  obtained,  processing  can  take  place  by  means  of  a 
computer  or  by  actual  optical  systems.  For  two-step  (holographic) 
processing,  optical  means  offer  a  natural  and  most  effective  approach. 

Even  with  one- step  processing,  it  seems  that  forming  an  image  might  be 
carried  out  by  optical  processing  and  display. 

llie  use  of  optical  experimentation,  both  in  simulating  target  return 
and  processing  to  reproduce  an  image,  requires  establishment  of  an  experi¬ 
mental  facility,  comprising,  at  least,  an  optical  bench,  an  associated  laser, 
and  photographic  equipment.  Also,  an  e7q}erienced  experimentalist  should 
be  obtained  to  operate  the  equipment. 

A  survey  was  made  of  appropriate  benches  with  associated  mounts, 
lenses,  spatial  filters,  alignment  sets,  together  with  gas  lasers  (which, 
in  fact,  are  quite  adequate  for  the  present  purposes).  Such  a  system  with 
sufficient  versatility  and  accuracy  would  cost  about  $10, 000. 00. 

In  addition,  microdensitometer  readouts  for  a  scanned  analysis 
of  the  hologram  or  image  patterns  can  be  obtained  on  a  rental  basis  for  about 
$40.00  per  run. 
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'  Considaring  the  importance  and  usefulness  of  the  experimental  aspects 

of  this  study,  the  adequacy  and  convenience  of  the  optical  mode,  and  the 
particular  low  cost  of  implementing  an  optical  set  up,  it  is  strongly  recom¬ 
mended  that  this  activity  be  initiated  and  supported. 
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SECTION  vn 
CONCLUSIONS 


MAIN  RESULTS 

The  puipose  of  the  present  studty  is  to  assess  the  feasibility  of  using 
large  microwave  apertures  for  obtaining  real-time,  three-dimensional  target 
images. 

Two  basic  approaches  are  considered:  (1)  one-step  processing,  which 
treats  the  array  as  a  lens,  and,  (2)  two-step  processing,  in  which  a  holo¬ 
gram  is  constructed  from  the  microwave  data  prior  to  forming  an  image. 

This  stu^y,  whose  preliminary  conclusions  are  summarized  below,  provides 
information  on  feasibility,  with  detailed  analytical  support. 

In  using  two-step  processing,  a  separate  reference  signal  is  required. 
This  can  be  either  synthesized  or  obtained  from  an  actual  reference  source. 
Very  accurate  knowledge  of  array  station  location  and  a  fully  coherent  array 
are  required  for  both  one-step  and  two-step  processing,  except  in  the  latter 
case  ndien  the  target  Itself  can  provide  the  reference  signal.  However, 
establishment  of  coherency  by  microwave  links  and/or  atomic  clocks  and 
determination  of  station  location  by  laser  ranging  appear  possible. 

When  a  station  spacing  based  on  the  spatial  frequency  content  of  the 
field  from  the  target  only  is  used  in  the  two-step  processing,  some  fold-over 
of  replicas  of  images  from  terms  of  the  hologram  on  the  desired  virtual  image 
of  the  target  is  caused.  The  fold-over  effect  can  be  dealt  with  without  altering 
the  station  spacing  by  defocusing  the  undesired  images.  Two-step  process 
would  be  dealing  with  at  least  a  hundred  stations  located  over  an  array  surface 
of  about  a  thousand  square  kilometers. 
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A  variety  of  approaches  are  available  to  reproduce  an  image  from  the 
sampled  microivave  data.  Three  approaches  have  been  noted,  with  particular 
emphasis  being  given  to  image  reconstruction  from  a  sampled  scaled  optical 
hologram.  A  calculation  of  the  effects  of  thermal  noise  is  given  for  tile  , 
sampled  holographic  process  and  SNR  is  determined  in  terms  of  scatterer 
properties  and  microwave  system  parameters. 

Althou^  the  lateral  resolution  can  be  under  a  meter  in  both  processes, 
the  longitudinal  resolution  is  much  poorer,  by  approximately  a  factory  of 
forty,  due  to  the  fact  that  the  range  to  the  target  is  at  least  about  ten  times 
the  maximum  linear  extent  of  the  array. 

To  date,  approaches  considered  for  using  the  target  as  a  microwave 
reference  source  appear  impractical.  However,  use  of  the  target  as  a 
reference  would  be  advantageous,  because  it  would  remove  the  strict  require- 
ment  on  accuracy  of  the  knowledge  of  station  location.  Such  accuracy  require¬ 
ments,  together  with  the  large  minimum  number  of  stations  and  the  poor 
longitudinal  resolution,  wei^  against  the  adoption  of  the  two-step  holographic 
method  without  the  use  of  the  target  as  a  reference  source. 

In  either  method  the  effects  of  off-zenith  target  position  do  not  ^[>pear 
to  be  a  serious  problem. 

Inasmuch  as  the  two-step  processing,  requiring  a  filled  array,  does 
not  yield  a  good  three-dimensional  virtual  image  because  of  the  poor  longi- 
tundinal  resolution,  it  may  be  advisable  to  explore  relaxing  the  three- 
dimensional  image  requirement,  thus  permitting  the  use  of  thinner  arrays, 
such  as  a  Mills  cross,  for  the  formation  of  a  two-dimensional  image 
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possessing  similar  transverse  resolution  capability  as  the  filled  array  con¬ 
sidered  for  use  in  the  two-step  processiiig. 

Application  of  such  arrays  together  with  possible  greater  diffusivity 
of  illumination  can  produce  the  advantages  of  improved  image  quality 
and  real-time  viewing  at  microwave  frequencies.  To  fully  exploxe  such 
advantages  in  image  quality  and  usefulness,  the  need  for  optical  range  data 
remains  an  Important  requirement. 

FUTURE  EFFORTS 

Now  that  the  various  major  advantages  and  disadvantages  are  under¬ 
stood,  analysis  and  experimentation  in  future  work  can  concentrate  on  partic¬ 
ular  areas.  Optimum  array  layouts  to  reduce  the  number  of  stations  and 
maintain  the  required  resolutions  will  be  considered.  The  enhancement  of 
image  quality  from  Illumination  diffusivity  needs  further  work.  Direct 
optical  methods  for  constructing  an  image  using  the  one- step  process  will 
also  be  considered.  Details  of  carrying  out  methods  of  obtaining  accurate 
station  location  determination  and  a  fully  coherent  array  can  be  established. 

In  addition,  methods  for  obtaining  a  reference  source  from  the  target 
will  be  investigated,  since  having  an  actual  reference  relaxes  the  coherence 
and  station  location  requirements. 

The  decision  as  to  the  superiority  of  the  one-or  the  two-step  process 
can  only  be  reached  after  both  approaches  have  been  fully  explored. 

Recognizability  of  the  image  in  various  target  environments,  together 
with  the  use  of  focusing,  reference  placement,  and  spatial  filtering,  need 
further  study*  The  formation  of  a  simulated  image  using  range  data  is  of 
prime  importance  here,  and  detailed  specifications  on  range  design  and 
experimentation  will  proceed. 
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APPENDIX  1 
CONTINUOUS  HOLOGRAPHY 

The  purpose  of  the  following  analysis  is  severalfold.  First,  it  is 
intended  to  provide  an  example  of  the  conventional  scalar  wave  equation 
treatment  of  the  formation  of  a  hologram  and  of  image  reproduction  there* 
from.  Second,  it  is  intended  to  provide  a  framework  for  the  study  of 
tolerances  that  must  be  observed  in  order  to  keep  within  acceptable  limits 
the  effects  on  the  reproduced  image  of  such  factors  as  a  rough  hologram 
recording  surface,  distortion  of  the  hologram  between  recording  and  repro¬ 
duction,  and  phase  errors  in  the  reference  signal  used  in  forming  the  holo¬ 
gram.  In  addition,  the  analysis  furnishes  the  basic  relations  that  are  re¬ 
quired  for  the  study  of  the  effect  of  the  location  of  the  source  of  the  reference 
signal  on  positions  of  the  reproduced  Images  and  on  depth  of  focus,  for  the 
study  of  the  magnification  that  may  be  achieved  by  the  system  as  a  function 
of  various  system  parameters,  and  for  the  study  of  the  reproduction  of  an 
image  from  a  "discrete"  hologram  synthesized  from  sampled  values  of  a 
continuous  hologram.  No  attempt  is  made  in  the  following  to  improve  on 
the  conventional  scalar  wave  equation  treatment,  except  for  the  simple 
provision  for  the  inclusion  of  various  error  terms  that  are  usually  omitted. 

FORMATION  OF  THE  HOLOGRAM 

In  Figure  1,  let  S  be  an  infinite  plane  tangent  to  the  earth's  surface 
H 

at  some  point  P  .  Let  a  coordinate  system  (x. ,  x  )  be  established  on  S 

1  1  ^  fi  ' 

with  P  as  origin,  and  define  the  unit  vectors  i.  and  i  along  the 

X  X  ^ 

X,-  and  X  -  axes,  respectively.  Erect  an  x  axis  perpendicular  to  S 
12  3  H 

with  the  direction  of  positive  x  being  to  the  left  in  the  figure.  At  a  dis- 

3 

tance  t  from  S  consider  a  plane  S  having  defined  on  it  coordinates  ({.,(  ), 

Xi  O  X  2 
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Figure  1 .  Hologram  Recording  Geometry 


witit  the  origin  of  coordinates  at  the  point  at  which  the  previously  defined 
X-  •>  axis  pierces  S  .  Ihe  -  and  -  axes  are  parallel  to  the  x  - 

9  0X2  X 

and  X  -  axes,  re«>ectlvely.  The  unit  vector  1  Is  normal  to  S  and  directed 

2  on 

to  the  left  In  the  figure.  On  the  plane  ,  let  there  be  defined  a  complex 
Amotion  D(i)  ,  which  will  be  taken  as  the  value  on  of  a  disturbance  propa¬ 
gating  in  space  In  accordance  with  the  scalar  wave  equation.  Harmonic  time 
dependeixe  of  the  form  e  Is  assumed. 

Define  the  following  vectors  for  arbitrary  x^  ,  x^  ,  x^  ,  ,  and  • 
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X  =  A,  X,  t- 


Xt. 


r  =  X  -t-  Ax  Xj 
f  =■  -?.  I,  ^  ^ 

?=  f  -  ?' 


It  is  assumed  that  the  hologram  is  to  be  recorded  on  a  surfkce 

Xj  =  J, 

The  reference  wave  used  in  forming  the  hologram  is  assumed  to  be  a 

spherical  wave  emitted  from  a  point  source  at  the  point  P_  located  at 

s 

*  ^3*3H  ■  '*“* 

s  ^  '' 

[10] 

If  a  quasi-monoc'jromatic  propagating  disturbance  of  wave  number  Ic,  takes 

on  the  values  D(0  on  S  ,  the  value  of  the  disturbance  on  the  surface  x  (i) 

Q  3'  ^ 

may  be  approximaved  bv[^^] 


where  r 


r|  .  If  it  is  assumed  that 


(1) 
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and 


over  the  region  of  8^  for  which  D(|)^0 ,  and  over  some  region  of  Interest 
of  Sjj ,  Equation  (1)  may  be  simplified  to 


>1 


Now  Introduce  an  approximation  for  r. 


(2) 


U  all  the  terms  except  the  first  in  the  above  bracket  are  small  conq)ared  to 
unity,  the  binomial  e]q;>anslon  theorem  yields  the  following  approximation 
for  r.  ^ 

....  V  (3) 


Assume  the  reference  signal  received  at  r’  is 


=■ 


t 


(4) 


where  B  is  a  constant.  Making  that  approximation  for  r^  that  corres¬ 
ponds  to  the  approximation  of  Equation  (3)  for  r  , 
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I  r„  /, »  _;i_ ) 

■LlJ-h.)  ' 


The  total  disturbance  ^^(r*)  reaching  the  point  r'  is  tL\e  sum  of  that 
due  to  the  distribution  D(i)  on  and  the  reference  signal.  Thus 


Now  in  the  phase  factors  in  Ekjuations  (2)  and  (4),  substitute  the  approxi¬ 
mations  for  r  and  r  ;  however,  in  the  denominators  of  those  equations, 

it 

which  are  slowly  varying,  substitute  the  more  crude  approximations  I  and 
I  +  ^  •  respectively.  Also,  in  the  phase  factor  of  the 

reference  signal,  incorporate  a  term  ik6  ,  where  6  m  6^)  ,  which  may  be 
used  later  to  represent  an  error  in  the  synthesis  of  an  artificial  reference 
signal.  Substituting  the  resulting  expressions  for  ^^(r')  into  Equation  (6) 
there  is  obtained 

f  -  Afj  ^  .'IC ..  LI  ,  ifi\  X,  ,1 


ti  er  r  r*  if  A  . 


The  power  flux  incident  at  r*  is  proportional  to 


!('')  = 'Tv (7 


89 


M68-2 


in  which  the  asterisk  Indicates  the  coii4>lex  conjugate.  Substituting  Equation  (7) 
into  Equation  (8)  yields 


I’ltTe  '  ^  >/'  / 


olijolf  J 

tit  ^  jIj}  ^ 


) 


JrJf' 


Air*  -i  ^  fw -<  *-  -XfAjUll.  L 

'  Uff}*  '  *  t-t  7i‘  /<(f 

'Vi.  r.A  -A,  ,_i!l  ,  ^.1^-7, n 

7 *-  e  J 

»  j  /i<  r  .1  ' 


lit  - 


[oVfi. 


%  ^ 


lllLil) 

/  .- 


8 

Cj  f  r,J'- 

o: 

The  quantity  !(?')  is  proportional  to  the  power  flux  density  on  the 
hologram  plane.  In  the  optical  case,  if  this  power  flux  density  distribution 
had  fkllen  on  a  piece  of  film,  the  film  could  be  processed  so  that  the  film 
density  would  be  proportional  to  the  logarithm  of  the  flux  density.  The 
result  would  be  an  optical  hologram. 

Having  obtained  the  hologram  corresponding  to  Equation  (9),  the 
task  is  to  reconstruct  an  image  of  the  original  "object"  -  in  the  present 
case  to  reconstruct  the  disturbance  D(|) .  For  ease  in  discussing  the 
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reooBBtruotioii,  let  Equation  <9)  be  written 


1(7')  ^  j, 


in  which  through  repreaent  the  four  terms  of  Equation  (9).  Ihst  ia, 
mairing  obvious  simplificationa,  and  using  r'  ■  r'  to  permit  replacing  f' 
in  the  argument  of  I  by  x  . 


fjrVV 


D(f)D*‘(r)^ 


.;A  Lillilj  illlzJll' 

A  "lA 


'if  Jij 


- iii _  f""  ‘  U{jfTy<i)  J^Tjix 

^irX(Jffr^) 

•  \o(s)e  ^  i-''  At 


{zrTrJrl 


vrj(/,  i,k) 


-  *  -illC  4 


o(f)e  (IS) 

> 


B  B* 

(j  *■  Tj,j^ 
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RECONSTRUCTION  OF  THE  IMAGE 

An  Image  may  be  reoonatructeo  by  passing  coherent  14^  through  the 
hologram.  For  generality,  It  will  be  assumed  that  a  spherical  coherent 
wavefront  is  used  for  reconstruction,  as  shown  in  Figure  2. 


Figure  2.  Reconstruction  Geometry 


Let  the  illuminating  radiation  be  of  wavenumber  k*.  The  hologram  will  be 
assumed  scaled  down  in  size,  so  that  one  unit  of  distance  in  the  S„ .  plane 
corresponds  to  m  units  in  the  plane.  The  point  at  x  -  0  in  the 
plane  will  be  assumed  to  lie  at  u  »  -u .  in  the  S  plane,  and  it  will  be 
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asBumed  that  the  hologram  has  beet'  distorted  somewhat,  so  that  the  point 

that  was  at  x  on  the  S  plane  maps  into  the  point 

H 


Ms  M  (J)  S 


(16) 


in  the  S.. .  plane, 
nd 


Here  A  =  A(u)  ,  and  A(-u .)  =  0  . 

d 


As  a  reference,  take  as  zero  the  phase  of  a  ray  that  has  Just  passed 

througdi  the  center  of  the  lens  in  Figure  2.  Tite  phase  of  that  ray  at  the 

focal  point,  P  ,  of  the  lens  will  then  be  k'f,  where  f  is  the  focal  length  of 

the  lens.  Now  consider  that  ray  through  the  focal  point  that  pierces  the 

hologram  at  point  P  located  at  u  .  Then  the  distance  from  P  to  P„  is 
H  Xf 


[> 


V 

I 


1 5  I' 


I 


y 


X 

V 


or,  retaining  the  quadratic  term  in  |u|  in  a  binomial  expansion, 


I  M 
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The  phase  of  the  disturbance  propagating  along  the  abovementicmed  ray  is 
then,  at  the  hologram  surface. 


Now  let  r^  be  a  vector  from  a  point  on  the  plane  to  the  point  .  Then 


''v’*  =  1 7;  |\  /  5  j  '*■  f  ^ 

“  I '  ^  —  I 


V 

V 


and,  retaining  quadratic  terms  in  |  or  -  u  |  in  a  binomial  eiqpansion, 


Using  Equation  (2),  the  disturbance  at  the  point  a  on  the  plane  may  be 
calculated  as  an  integral  over  the  hologram.  At  the  rij^  hand  side  of  the 
hologram,  the  disturbance  ^  is 


ACf 


r 


^  C<; 


where  A  is  the  constant  of  proportionality  relating  the  power  flux  density  I 
to  die  fraction  of  the  incldmit  field  intensity'*'  transmitted  toe  hologram  in 

*  hi  this  work,  intensity  refers  to  electric  field  strength. 
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the  reconstruction  process,  Ct/z^  is  the  intensity  of  the  coherent  tiinmiTiatinn 
incident  on  the  hologram,  and  C  is  the  intensity  of  the  plane  ivave  incident  on 
the  lens.  In  the  MKS  system,  the  units  of  A  are  meter  /volts  .  Now,  the 
disturbance  at  or  is,  using  Equation  2, 


-  d  f 


Setting  r^  »  in  the  denominator  of  the  above,  and  using  the  usual 
binomial  expansion  approximation  for  r  in  the  phase  factor,  there  results 


^  A  ( 

XTT  ♦  ,  ^  ^ 


♦  -f 
e 


) 


'' 


V 


X 


// 


<Jir, 


(16) 


Equation  (16)  is  linear  in  the  terms  throu^  comprising  I .  Let 
us  determine  the  contribution  to  ^^(a)  of  each  of  these  four  terms. 
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(a)  Contribution  from  .  Let 


■A-i^ 


Kft  w  1 


(17) 


Now  aasume  that  flie  region  of  4  in  the  plane  S^  over  which  D({)'/  0 
is  suttioiently  restricted  so  that  the  plane  S  lies  in  the  Fraunhofer  region 
of  aqy  radiation  arriving  at  from  the  distribution  D(l) .  That  is,  for 
all  i  for  which  D(i)#  0 , 


klJt 


<<  I 


(18) 


Further,  assume  \^X^X  >  *• 

The  restriction  that  die  hologram  plane,  8  ,  lie  in  die  Fraunhofer 

H 

region  of  the  scattering  complex  is  not  essential.  If  the  condition  of  Equa- 
tion  (18)  is  not  met,  the  phase  factor  eaq>(-ik  |  { |  /2I)  may  be  groiqped  with 
the  disturbance  D(i) ,  and  the  entire  analysis  carried  throuj^.  The  recon¬ 
structed  image  will  then  contain  this  same  phase  foctor,  but,  in  optical 
viewing,  it  is  the  squared  amplitude  of  the  field  that  is  important,  and  this 
squared  amplitude  is  unchanged  fay  the  phase  factor.  The  entire  resulting 
procedure  is  termed  Fresnel  holognqilqr. 

Define 


96 


M68-2 


The  plane  ie  preeumed  infinite,  altboui^  D(4)  ia  non-sero  over  (Moly  a 
finite  region  of  .  Letting 

^  *  — I - c  / 


and  uaii^  the  Fraunhofer  c<mdition,  Equation  (18).  Equation  (17)  beconuw 


C  ‘  ^  ^  ^  'tf 


Now  impose  the  focusii^  condition 


Equation  (21)  becomes 


.  (22) 


Assuming  A  <=  0  (this  is  convenient  but  not  necessax]^)  the  above 
integral  will  be  expressed  as  a  convolution.  Let 


L  u  t  uj 

w  =  t  - - i- 


97 


1168-2 


Vi)’  Hf'  Je"  S  )i- 


-  fi 


r 

b’(Z)  t>  (7- 


(23) 


(b)  Contribution  from  I  .  Let 

A 


j.  oy  -SI^Fjn 


e 


^"Ti;  ^Mtr„^r  ...  r 


_  -i4|r.-l 

0/lje  'f 


A 


ifi'' 


•  .  *1  /A -f/' 

.  .  ^-<  "  “TIT- 

^  T  ^  ^ 


(24) 


Now,  using  the  Fraunhcrfer  condition 


iilC 

V  jf 


<  <  i 
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:^*.i  .acrt 


r.s  *  *  \J. 
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and  assuming  that 


ilkjl  lx  - 


(26) 


and 


ilrill  k-r^ 


(26) 


which  impose  requirements  on  ,  Equation  (24)  becomes 


HTTTTrZT' 


<  ^J-i)  fn 

^  *■  ftK 


xn  P  /  km« 


[i. 


- 


i-)j- 


(27) 


where  is  the  factor  multiplying  the  double  integral  in  Equation  (24). 
Now  note  that 


llA  t  UJ-  Ij'' 


(28) 
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To  obtain  the  focusing  condition  for  $  ,  In  the  exponent  equate  to  zero 

«  »2 
I  1^ 

coefficient  of  i 


Using  Equations  (28)  and  (29),  Equation  (27)  becomes 


..As 


X  e 


*  e 


'  ■  -  ^  •  (m;  -aJ  ~  /. 


In  the  case  in  \i«diich  the  errors  6  ,  and  A  ,  are  zero,  this  is 
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Lettins 


•  xtFx  f 


,  this  becomes 

^  ;a  L 

♦  -"l.V 


ju*-  t,^)  Js(^)j:^^ 


so 


.4 


..  <«^ir 

XiOfTit) 


.1 


X 


■A*"'  r,^  \7 

j]  (30) 


Representing  the  entire  coeffici«it  of  D  by  F  ',  the  esqpression  for  ^  (^ 

ta  'ss 
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Equation  (31)  Indicates  that  at  the  value  of  satisfying  Equation  (29). 
there  will  appear  an  image  of  the  original  distribution  on  .  The  image  is 
magnified  and  translated,  and  may  be  a  real  or  virtual  inu^e  depending  on 
whether  is  greater  than  or  less  than  zero,  respectively. 

(c)  Contribution  from  Ig  .  Let 

**  i  ^  ^  ^  I )  L  ^  ^  ^  If  i  V  t  t  Ti  r 

.  e  ~  *J  ■'  ~  -rjTr^fr^ 

*e  TTTTTTT  ^  ^ 

J  4t4:z 


and  let 


Using  the  Fraunhofer  condition. 


<  <  I 


i 


M68^ 


and  amiunlng  fiviatloBS  <25)  and  <26)  are  valid,  Eqttttton  (33) 


c  O'" 


^  4  A  M*  * 

xf  -/♦  T,„ 


VTT  D 


Using  Equation  (28),  insertinf  tiie  focusing  condition 


*ii-  •  tj  -  ,  0 

'  '  ' 


there  results 


777771^  " 


I 
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Again  conalder  the  case  in  which  the  errors  6  and  A  ,  are  zero. 


Then 


^  («  I  g  I  / 


Letting 


,e  ■<(/*'..)  f  Y<- 


W  *.  g  ■  ■  I  -  ^ 


-  ■^■>  ;  fi  M4'  ?rt 

^t,,  s  ^.n-  f ,  e  Tt*)  x^T,a 

J>  r  -(r^-'-t)  (^ '  . <i»,'.5.'i!L ) 

'^|e  '  '  '  -'•I* 

»  ^‘>^J 

.Jk^r  iikl/k  iiiiii.  f  \ 

^  I  Pj  f  >«-/>ii  I  >•  AU*r,A)f 


.  p-  /  ,  >(  F  “  -. 

>  y  I  — :  *•'  + - -  j.  '‘' 

^ /  Tr^ 

<  -v  '*  *  *r  p  / 
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Appljring  this  conditim, 


That  Is, 


li 

JT¥ 


SO  that  in  flie  z-plane  In  vdiich  z  satisfies  the  focusing  condition,  (o) 
is  a  delta-function. 

In  the  foregoing,  the  plane  ,  on  (or  near)  which  the  hologram  is 
formed  has  beat  assumed  to  be  infinite  in  extent.  This  necessitates  the 
violation  of  the  assumption  made  in  going  from  Equation  (1)  to  Equation  (2), 
that 


r 

In  praotloe,  of  course,  the  hologram  is  of  finite  extent.  The  finite  size 
of  ttie  hologram  may  be  taken  into  account  in  the  foregoing  mathemattos  by 
multiplying  the  ri|^t  side  of  Equation  (9),  and  hence  the  rii^t  sides  of 
Equations  (11),  (12),  (13),  and  (14),  by 
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R..*  «,x  )  R«<4  (7-)^ 

where  the  function  Rect  (x)  is  defined  to  have  the  value  unity  for  -l<x  $  1 » 

Mrp  olherwl«.  «id  Xj  «id  ar.  conataat..  Tbl.  lactor  «ffl  rlda 

through  the  entire  analysis,  with  suitable  substitution  for  and  x^  in 

terms  of  u^^  and  u^  ,  and,  in  the  end,  expressions  for  ,  and 

Ip,  will  be  obtained  that  are  convolutions  of  the  expressions  that  have  been 
^84 

obtained  above  for  these  quantities  with  the  Fourier  transforms  of  the  multipty- 
ing  Rect  functions.  The  Fourier  transform  of  the  Rect  function  is  a  sinx/x 
function  diat  approaches  a  delta-function  as  and  tend  to  infintty. 

Thus  for  large  but  finite  values  of  X^  and  X^  ,  the  effect  of  the  finite  size 
of  the  hologram  is  to  smear  the  Images  that  have  been  predicted  in  the  above 
treatment,  with  the  resulting  resolution  being  determined  in  the  usual  way  fay 
the  ratio  of  and  X^  to  the  operating  wavelength. 

SOME  PARAMETERS  DESCRIBING  THE  HOLOGRAPHIC  PROCESS 
Magnification 

From  either  Equation  (31)  or  Equation  (37)  describing  the  images 

corresponding  to  the  quantities  I  and  I  from  Equations  (12)  and  (13), 

z  3 

respectively,  it  is  ai^>arent  that  in  the  case  that  a  real  image  is  formed, 
a  unit  length  in  the  object  plane  S^  corresponds  to  a  length  mkz^/llc'  in 
the  image  plane  S^  .  Accordingly,  the  magnificati<m  of  the  holognq)hically 
produced  image  will  be  defined  as 
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►VI,,  X 


a. 

A  A' 


(40) 


The  eubecript  ()  on  m||  stands  for  parallel,  as  a  reminder  of  the  fact 

that  the  magnification  corresponding  to  directions  lying  in  the  plane  S 

0 

is  different  from  thnt  corresponding  to  the  direction  perpendicular  to  that 
plane.  To  obtain  the  magnification,  n^  ,  corresponding  to  the  dlrectior 
perpendicular  to  K  ,  consider,  for  example,  the  focusing  condition  C  : 


♦  r;  . . '  «  0 


The  magnification  mi  is  found  by  inquiring  what  change  in  z  is  required 

to  bring  into  focus  an  object  at  I  dl  ,  assuming  the  hei|d*t  above  8 

H 

of  the  reference,  I  ,  to  remain  constant  Letting  I  ^  t  +  , 

3R  A  3R 

O  may  be  rewritten  as 


(;  = 


(41) 


Setting  do  ■  0  for  simultaneous  variation  of  z  and  I  ,  there  results 


ltt8-2 


Hence, 


and 


"X" 


(42) 


(43) 


Resolution 

The  resolution  capability  of  a  holographic  system  depends  on  whether 
the  two  points  to  be  resolved  lie  in  a  plane  parallel  to  the  array  aperture  or 
on  a  line  perpendicular  to  that  aperture.  These  two  resolution  capabilities 
will  be  termed  parallel  and  perpendicular  resolution,  respectively. 


Parallel  Resolution 


As  described  near  the  end  of  the  preceding  section  of  this  Appendix, 
the  use  of  a  holographic  recording  plane  of  finite  extent  leads  to  expressions 
for  images  that  are  convolutions  of  the  e^qpressions  that  have  been  obtained 
for  an  infinite  aperture  with  the  Fourier  transform  of  the  actual  finite  aper¬ 
ture  shape.  In  the  case  of  a  square  aperture,  the  function  with  which  the  true 
image  is  convolved  is  a  two-dimensional  sinx/x  function.  Consider  the 

result  for  the  case  of  a  point  source  in  the  object  plane;  i.  e. ,  D(i)  is  a 

_  -1 

6-function.  By  Equation  (17)  ,  setting  D({)  -  (2r)  , 


A  ^  T,fi  J 

J 


(44) 
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whare  the  aperture  in  the  hologram  plane,  8^  ,  is  t^k^n  to  be  a  square  of 
length  2L  on  a  side.  Let 


V  = 


fit 

J  f 


Carrying  out  the  Integration, 


(45) 


in  which  and  are  the  components  of  v  along  the  and  axes, 

respectively.  If  two  point  sources  are  present,  they  will  be  said  to  be 

resolved  when  the  first  zero  of  the  hmction  p-  (c^  for  one  source  occurs  at 

182 

the  peak  of  the  function  for  the  other  source.  If  the  sources  are  dis¬ 

placed  along  the  direction,  this  corresponds  to  a  variation  in  the  v^ 

variable  of  v, .  such  that 
Id 


*vt 


TT 

T  yv1 

1. 


nils  corresponds  to  a  variation  in  cc^  of  ,  where 


= 


ik'L  ^ 
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which,  in  turn,  corresponds  to  a  displacement  in  the  object  plane. 
Using  Equation  (40) 

I  -  ^ 


(Strictly  ^'peaking,  this  result  would  be  modified  by  a  fkctor  (A  <^2  ii  the 
displacement  of  the  two  point  sources  had  been  in  any  of  the  directions 
^(Ui  ^  u^/  .  This  change  will  be  ignored  ,  inasmuch  as  the  definition  of  resolu¬ 
tion  is  somewhat  arbitrary  anyway. ) 

Perpendicular  Resolution 

The  perpendicular  resolution  can  be  determined  by  inspection  of 
Equation  (27)  or  Equation  (34).  Consider  Equation  (27)  and  assume 

u .  -  0  and  A  =  0  -  0  .  Assume  that  a  point  source  lies  at  i  -  0 

®  -  -1 
in  the  object  plane,  .  Then  D  (kmu/l)  =  (2ir)  ,  and 


SI 


(47) 


For  G  =  0  ,  the  behavior  of  il>.  has  been  investigated  in  the  above  deter- 

^2 

mination  of  the  parallel  resolution.  Now  assume  that  v  >  0  and  investigate 
\l>  as  a  function  of  z  .  That  is,  now  investigate  the  (possibly  virtual) 

o2  ^ 

reconstructed  field  along  an  axis  perpendicular  to  the  focal  plane  of  the 
image  and  passing  throue^  the  center  of  the  image  in  that  plane. 


Ill 
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Equatton  (47)  becomes 


t  u. 


0 


(48) 


Note  that  if  G  =  0  , 


IL  /r.  I  U 


Equation  (48)  may  be  rewritten  as 


where  the  (Fresnel)  Integral  In  the  above  equation  is  tabulated  .  Taking 

the  magnitude  of  ^ 

%2 


I  IV  I 


/C  ~  < 


Vm) 


S‘ 


(49) 
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The  functions  C  and  S  are  tabulated  in  Reference  12.  The  function  of  G, 
U  P  ,  has  no  zeroes.  At  G  =  0,  P  has  the  value  1^,  I  ^  .  As  G 

I  ^2l  I  %2l  1  IS2|o 

becomes  larger,  the  quantity  in  brackets  in  Equation  (49)  grows  to  a  value  of 
the  order  of  unity,  then  decreases  to  a  value  of  the  order  of  0.05,  after 
which  it  increases  again  and  oscillates  with  decreasing  amplitude  about  its 
asymptotic  value  of  0.25.  As  a  resolution  criterion,  it  will  be  considered 
that  a  second  point  source  will  produce  an  image  that  can  be  resolved  firom 
that  of  the  first,  if  the  center  of  the  second  image  falls  at  a  value  of  z 


such  that  ip. 


is  0.1.  That  is. 


0,  I  = 


TT'"  ^ 

7^ 


This  condition  is  satisfied  when 


f  T 


^  I.  -T 


From  Equation  (41) 


c «  — 


80 


ISVvl 


"Pi' 


This  displacement  in  image  space  corresponds  to  a  displacement  In 
object  space,  where 


I 


(50) 


Thus  the  perpendicular  resolution  is  usually  much  poorer  than  the  parallel 
resolution. 
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APPENDIX  2 

WAVEFRONT  RECONSTRUCTION  FROM  SAMPLED  DATA 

In  performing  microwave  imaging,  whether  by  a  holographic  technique 
or  by  a  method  analogous  to  optical  image  formation  with  a  lens,  the  scattered 
field  is  distributed  over  such  a  large  region  that  it  is  generally  impractical  to 
record  the  field  on  a  continuum  basis.  Consequently,  it  is  necessary  to  resort 
to  sampling  the  scattered  field  and  to  some  method  of  reproducing  an  image 
from  such  sampled  data.  If,  in  the  case  of  holographic  imaging,  an  actual 
reference  signal  were  present,  it  would  be  desirable  to  measure  directly  the 
amplitude  of  the  resultant  of  the  reference  field  and  the  field  scattered  from  the 
target,  because  by  so  doing  the  requirement  on  knowledge  of  station  location 
would  be  minimized.  However,  because  no  actual  reference  signal  of  con¬ 
trolled  location  is  available  in  any  of  the  microwave  imaging  techniques  that 
have  yet  been  devised  for  use  in  this  investigation,  all  methods  for  micro¬ 
wave  imaging  considered  as  yet  in  this  study  assume  that  the  amplitude  and 
phase  of  the  scattered  signal  from  the  target  are  recorded  at  each  station. 

If,  then,  it  is  desired  to  reconstruct  the  wavefront  at  optical  frequencies  from 
a  hologram,  the  amplitude  and  phase  that  a  reference  signal  from  a  source  at 
a  given  location  would  have  at  each  of  the  observing  stations  can  be  computed 
(assuming  station  location  is  known  with  sufficient  accuracy — see  Section  IV), 
and  the  computed  reference  signal  can  be  combined  with  the  measured 
scattered  signal  to  find  the  resultant  field  amplitude  that  would  have  been 
observed  at  each  station,  if  an  actual  reference  signal  had  been  present.  This 
measured  amplitude  (or,  more  precisely,  the  square  of  this  amplitude) 
furnishes  samples  corresponding  to  those  of  a  hologram  measured  at  the 
various  stations. 
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In  the  two  techniques  of  wavefront  reconstruction  described  in  the 
following  under  Methods  1  and  2,  the  abovementioned  hologram  samples  are 
used  to  form  a  "discrete"  hologram,  scaled  down  suitably  for  optical  repro¬ 
duction  (see  Example)  that  consists  of  a  set  of  square  holes  in  an  opaque  screen, 
one  square  for  each  sampling  station,  with  each  square  having  constant  light 
transmission  over  its  area.  Methods  1  and  2  are  alternate  methods  of  recon¬ 
structing  a  wavefront  from  a  sampled  hologram,  assuming  the  initial  sampling 
is  done  at  intervals  over  the  earth's  surface  suitable  for  recording  data  from 
a  target  of  some  specified  maximum  size.  Of  these  two  methods,  the  first  is 
considered  superior  because  of  simplicity  and  reduced  sampling  requirements. 

In  addition  to  Methods  1  and  2,  the  following  presents  still  another 
wavefront  reconstruction  technique.  Method  3,  not  making  use  of  a  hologram 
at  all,  but  reconstructing  the  wavefront  directly  from  the  sampled  values  of 
scattered  field  amplitude  and  phase,  ’^e  requirement  on  the  spatial  interval 
at  which  the  scattered  field  is  sampled  is  the  same  in  Methods  1  and  3. 

Methods  1  and  2  require  that  the  hologram  be  illuminated  by  a  coherent  plane 
wave  and  that  the  amplitude  of  the  illumination  passing  through  each  square 
of  the  hologram  be  controlled  by  varying  the  transmission  of  the  squares. 

Method  3  requires  that  both  amplitude  and  relative  phase  be  controlled  of 
each  of  a  number  of  beams  (the  number  being  equal  to  the  number  of  sampling 
stations).  Methods  1  and  2,  of  course,  each  require  the  synthesis  of  a  refer¬ 
ence  signal,  which,  in  turn,  requires  a  knowledge  of  sampling  station  location. 
Method  3  does  not  make  use  of  a  reference  signal,  but  places  equally  stringent 
requirements  on  knowledge  of  station  location  in  order  to  determine  the 
appropriate  phaser  of  the  beams  used  in  the  direct  wavefront  reconstruction. 
Method  1  is  preferred  over  Method  3  as  well  as  over  Method  2,  because  it 
appears  to  be  the  simplest.  Because  Method  1  is  currently,  at  least,  the 
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preferred  technique  for  reconstruction,  it  will  be  presented  in  some  detail, 
while  Methods  2  and  3  will  be  dealt  with  more  briefly. 

METHOD  1;  DIRECT  RECONSTRUCTION  FROM  CHECKERBOARD 
HOLOGRAM 

In  the  following,  the  hologram  will  still  be  considered  of  infinite  extent 
(the  effect  of  using  a  finite  hologram  may  be  determined  by  the  procedure  de¬ 
scribed  near  the  end  of  the  preceding  Appendix),  but  it  will  be  assumed  that  the 
value  of  the  function  I  of  Equation  (9)  of  Appendix  1  is  known  only  at  the  points 
*lj  ’  *2n  square  grid.  The  points  x^^  and  x^^  map  into  the  poi  its 
u, .  ,  u„  in  the  S„ ,  plane  Iw  Equation  (15)  of  Appendix  1.  In  the  following, 

ij  £ti  llQ 

the  error  functions,  6  and  A  ,  will  be  taken  to  be  zero,  inasmuch  as 
it  is  the  effect  of  using  a  sampled  hologram  that  is  to  be  investigated.  Further, 
for  the  present  purpose,  only  one  of  the  terms,  I^  ,  of  the  four  terms  in 
Equation  (10)  of  Appendix  1  will  be  considered.  The  displacement  u^  will 
be  taken  to  be  zero,  for  simplicity. 

An  obvious  way  in  which  to  attempt  to  reconstruct  an  image  from  the 
sample  values  of  the  hologram  specified  above  is  to  make  a  hologram  having 
its  transmission  constant  over  square  cells  of  length  2b'  on  a  side,  where 
2b'  is  yet  to  be  specified,  with  the  cells  centered  in  the  scaled  hologram  at 
the  points  corresponding  to  the  location  of  the  receiving  stations  in  the  original 
hologram  plane.  The  constant  transmission  over  each  cell  is  that  value 
determined  by  the  measurement  at  the  corresponding  receiving  station.  Re¬ 
construction  of  an  image  from  such  a  discrete  hologram  will  now  be  con¬ 
sidered. 
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Making  use  of  Equation  (12)  of  Appendix  1,  setting  6  =  0,  the  second 
term  in  the  discrete  hologram  may  be  written 


where  x^^^  =  2jb  and  x^^  =  2nb  .  K  this  hologram  is  mapped  into  the 
u- plane,  using  Equation  (1 
hologram  in  the  u-plane  is 


u-plane,  using  Equation  (15)  of  Appendix  1,  with  u,  and  A  both  zero,  the 

d 


(1) 


Using  Equation  (16)  of  Appendix  1  for  the  reconstruction. 


'**4 


(2) 


Here  mu,,  =  x,,  and  mu  =  x„  .  Let 
Ij  Ij  2n  2n 


a  * 


I 


(3) 
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Consider  now  the  integral  over  u  .  The  following  integrations  are  involved 


It  will  now  be  assumed  that  the  distance  b'  is  sufficiently  small  so  that  for 

I  I  b*  I  I  b* 

all  u,  and  u„  such  that  u,-u,.  and  u„- u„ 

1  2  ‘  1  Ij'  m  '  2  2n'  m 


and 


a  1 


(4) 


(5) 


for  all  j  and  n  .  In  the  above  integrals,  the  quadratic  phase  factors  may 
then  be  considered  constant.  The  result  of  the  integrations  over  u^  and  u^ 
is  then 


.k'  -  V  - 


.  ^  ill: 


kV,  V) 

JIm-k  \  I 


k  <  f 


where 


t. 


J,'*' 


’» 


A 

U 


f  'A 
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Substituting  the  result  of  the  u  integration  into  Equation  (2)  , 


^  ^  -tt: 

'  TTa,*^  k.' oft  k>'  ^ 

tt.  ii 

’  b'  *  W'. « I*”  ;  k'  Of  • 

ZiK  J  4.*» '  -  »«.  — .  4i^  .  . 

•t  ^  ^  (6) 


Now  consider  only  the  summation  in  Equation  (6),  and  substitute  from  Equa¬ 
tion  (12)  of  Appendix  1,  with  6=0,  and  as  usual  assume 


liil" 

R  |x,))x /<1 

2  4*  ^ 


« 1. 


The  summation  becomes 


-i^SsK' 


.  ^  CiJJlizjlii*  4  V  V  .  «  L  ^  ^ 
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Now  let  it  be  assumed  that  D(4)  is  zero  except  possibly  over  the  square 
region  in  the  plane  , 

-Lo  ^  i  ^  b 

-  Lb 

in  this  case  D(|)  may  be  considered  to  be  repeated  periodically  over  the 
plane  ,  and,  in  the  square  interval  specified  above,  may  be  repre¬ 
sented  by  a  Fourier  series  as  follows, 

—  ‘  (t,j  -*  ii"** 


Here 


(9) 


’TT'Tt 

Ld 


) 


(10) 


and  the  last  step  follows  from  the  vanishing  of  D(|)  outside  of  the  square 
region.  Thus 


r(T)==  \ 


(11) 
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Now  insert  in  Equation  (7)  the  focusing  condition 


to  obtain 


.  kb 


i!hi" 


-f 


-tM 


Ji-n 


(  k'^ 

-  k^lR  N 

•  27fD 

(12) 


277  D  (13) 


Note  that 


s  r.  K'v.' 

■«  “1 —  — 

—  t,  ^  ^  ^  4-  ti. 

A  I 

and  let  b  be  such  that 

7L  -  b  (14) 

i-i>  ~  1 
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Then,  using  Equations  (11)  and  (14),  Equation  (13)  may  be  written 


or,  again  using  Equation  (14), 


^  i 


-  .Tm^n  UL  d  -  iii.  +  iit.  \ 


(16) 


Now  substitute  from  Equation  (16)  for  the  summation  in  Equation  (6). 


-  A  k  B  C  f 


aV, 


“ZTii 


«%!» 


*  «r.  k» 

^  a  k 


h*  w 

—  1. 


-  :  *  %j«  •  '  ^  JIs!— 
<e  ‘«’U 


,  \  — 

^  k 


(17) 


Thus,  if  sample  *  the  hologram  is  done  at  a  spacing  of  2b,  where 
b  is  not  larger  than  the  disumce  determined  Equation  (14),  and  if  b' 
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satisfies  the  inequalities  (Equations  4  and  5),  which  are  necessary  to  permit 
focusing,  Equation  (17)  indicates  that  an  image  of  the  distribution  D(i)  can 
be  produced  from  a  sampled  hologram. 

In  Equation  (17)  there  are  two  sinx/x  functions  multiplying  the  desired 
image  function.  To  minimize  degradation  by  this  intensity  modulation,  the 
desired  image  should  be  placed  between  zeroes  of  the  modulation.  This  can 
be  done  easily.  Another  source  of  image  degradation  lies  in  the  periodicity 
of  the  reconstruction  described  by  Equation  (1.3).  That  is,  the  image  de¬ 
scribed  by  Equation  (16)  is  repeated  periodically  in  the  a-plane,  with  period 
m||2Lj^,  where  m^i  is  the  previously  defined  parallel  magnification.  To 
prevent  serious  interference  with  the  desired  image  from  the  periodically  re¬ 
peated  reconstruction  fields  due  to  the  hologram  terms  ,  Ig  ,  and  ,  the 
system  should  be  operated  with  |4g^|  sufficiently  different  from  zero  so 
that  the  desired  reconstructed  image  (i.e. ,  that  due  to  either  I  or  I  ) 

(6  J 

can  be  isolated  by  focusing.  Because  of  the  focusing  condition  the  value  of 
that  is  used  sets  an  iqiper  limit  on  b'  through  Equations  4  and  5.  This 
u(^r  limit  is  not  a  source  of  difficulty  in  system  design,  however. 


METHOD  2 


In  this  method,  the  sampled  hologram  data  is  used  to  form  a  con¬ 
tinuous,  scaled  down,  rec<nistructlon  of  the  original  hologram,  which,  in 
turn,  is  used  to  reconstruct  an  image. 


Let  I  (u)  s  I(mu) ,  and  assume  I  (u)  to  be  known  at  the  points 
8  8 

(u,,  »U-_)  of  a  square  grid,  where  u,.  and  u^  have  been  previously 
IJ  ■  2n  ij  2n 

defined.  As  before,  the  function  I^)  is  that  given  by  Equation  (9)  of 

Ai^>endix  1.  It  is  assumed  here  (in  distinction  to  the  case  in  Method  1) 

that  the  spacing  of  the  points  at  which  I  (u)  is  given  is  such  that  I  (Q)  could 

8  8 
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be  reproduced  exactly  by  use  of  the  sampling  theorem,  if  the  square  grid 
were  of  infinite  extent.  This  assumption  introduces  additional  questions 
inasmuch  as  the  hologram  contains  a  quadratic  phase  term,  the  frequency 
corresponding  to  which  increases  indefinitely  as  the  size  of  the  hologram 
is  increased.  However,  it  may  be  assumed  that  the  hologram  is  continued 
beyond  the  region  in  which  it  is  known  in  some  manner  other  than  that  in 
which  it  actually  continues.  With  this  assumption,  there  arise  questions  of 
what  sort  of  continuation  (perh^s  periodic)  might  be  assumed  in  order  to 
improve  the  representation  of  the  actual  hologram  by  sampling  functions, 
artificially  extending  the  hologram  beyond  the  region  in  which  it  is  known. 
These  questions  will  not  be  considered  here,  however.  Instead,  it  is  simply 
assumed  that  a  function  is  given  over  the  entire  u-plane  that  coincides  with 
the  hologram  over  the  known  region  of  the  hologram,  and  that  is  reproduced 
exactly  by  the  sampling  theorem  (at  least  over  the  known  region  of  the  holo¬ 
gram).  An  estimate  of  the  required  sampling  interval  will  be  given  later. 

Let  the  "discrete"  hologram. 


^  1 1 

be  formed  from  the  mt  Mred  values,  I  (u  ,)  =  I(x  ,)  ,  with  x  = 

B  n,J  n,j  n,j 

mu  .  .  Let  the  FourU  ransform  of  the  central  square,  i.  e. ,  that  corres- 

j  _ 

ponding  to  u^  ^  be  taki  optically,  using  a  lens  "f  to  f. "  The  resulting 
dishirbance  in  the  back  foe.  i<.ae  of  the  lens  is 
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where  ^(x)  is  defined  to  be 


H...  !.») 


<  *  U. 

t  A«4. 


*1-  f  Uk(1 


I. .<,(*) 


j  J  —  -  V  -rr*  -=^ 

TA*^.,.)  - 2= - 

rf.  w*  K*  w 


Inserting  Equation  (21)  into  Equation  (19) , 


^i")  e  ii-,  Xjly.,.) - - 12.. 

t  Tr~.‘  w,  b-  h'  »’ 

<4  ^ 

If  the  same  processing  is  performed  for  the  square  with  center  at  u  ,  , 

J 

with  the  axis  of  the  lens  moved  to  the  location  u  ,  ,  the  disturbance  in 

“•J 

the  w-plane.will  be 


t-.i  (-J  = 


^  no  *  f  nin 


>  .  .lb 
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The  sum  of  such  disturbances  over  n  and  J  Is  then 


•  "  ^  k‘W‘  /  u  \  v'k.'/  J 

Let  the  focal  length,  f,  of  the  lenses  used  in  this  processing  be  such  that 

the  zeroes  of  the  sin  x/x  function  appearing  in  the  summation  fall  at  the 

sample  points  w  ,  (=  u  , )  ,  spaced  by  2b/m. 
n.j  n,j 

Thus  it  is  required  that 

^  s  TT 

A  ^  > 

T  ^ 


f  » 


where  =  2  tt  /k'  .  Using  this  value  of  f ,  the  summation  on  the  right  of 

Equation  (24)  is  I  (w)  ,  by  the  sampling  theorem.  The  total  disturbance  in 
s 

the  back  focal  plane  of  the  lens  system  is 

e  Xj  (26 

where  Equation  (25)  has  been  used  in  Equation  (24).  Thus,  except  for  a 
constant  phase  factor,  the  unsampled,  continuous  hologram  has  been 
recreated  in  the  w-plane.  This  hologram  may  now  be  utilized  in  recon¬ 
structing  an  image  of  the  disturbance  that  gave  rise  to  the  original  hologram 
l(x)  . 
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Sampling  Requirements 

As  stated  earlier,  it  is  the  object  of  the  present  reconstruction 
method  to  reproduce  first  from  the  sampled  hologram  a  continuous,  scaled 
down  reproduction  of  the  original  hologram,  from  which  the  actual  recon¬ 
struction  can  be  made.  Hence,  it  is  necessary  to  sample  the  original  holo¬ 
gram  at  a  rate  sufficient  to  reproduce  the  highest  spatial  frequencies 
occurring  therein.  Inspection  of  Equations  (11),  (12),  (13),  and  (14)  in 
Appendix  1  reveals  that  the  spatial  frequency  spectrum  of  the  original 
hologram  may  be  represented  schematically  by  the  following  diagram. 
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In  order  to  obtain  spatial  separation  of  the  desired  image  from  inter¬ 
ference  terms  in  the  reconstruction,  |  |  should  be  chosen  at  least  large 

enough  so  that  the  frequency  ranges  shown  in  Figure  1  do  not  overlap.  With 
chosen  as  small  as  possible  without  overlap  occicrring,  the  spectral 

content  of  the  hologram  is  seen  to  occupy  the  range  (-  * 

V  IT  ^  ■’T 

requiring  a  sampling  rate  of  iLkiJL  ,  or  a  sampling  interval  of 

TT  f 


Li\ 
(  ^ 


-  ±Ji 


(27) 


This  is  one-fourth  of  the  interval  required  Method  1,  and  results  in  six¬ 
teen  times  as  many  sampling  stations  as  that  method  This  is  not  to  say  that  it 
might  not  be  possible  to  form  a  reconstruc  ted  hologram  (not  image)  as  in 
Method  2,  using  thinner  sampling,  and  then  to  separate  the  desired  image  in 
the  reconstruction  by  use  of  focusing.  However,  this  has  not  yet  been 
investigated. 

METHOD  3:  DIRECT  WAVEFRONT  RECONSTRUCTION 

R.  W.  Roig  of  the  MITRE  Corporation  has  proposed  a  wavefront  recon¬ 
struction  technique  that  does  not  require  the  use  of  a  hologram  at  all.  A 
hologram,  of  course,  serves  as  one  method  of  preserving  phase  information 
concerning  an  incident  wavefront.  In  the  microwave  case,  both  amplitude 
and  phase  would  be  recorded  at  each  receiving  station  in  any  imaging  method 
not  using  an  actual  reference  beam.  The  phase  information  about  the  received 
wave  is  then  alreacty  recorded,  and  Roig  has  pointed  out  that  it  is  not  necessary 
to  use  a  hologram  at  all  in  order  to  reconstruct,  at  another  frequency,  a 
scaled  approximation  to  the  orig^inal  waveffont. 
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One  significant  disadvantage  of  the  Roig  reconstruction  method  is  that 
it  requires  not  only  control  of  the  amplitude  of  a  discrete  set  of  coherent 
sources  of  optical  illumination,  but  control  of  their  relative  phases  as  well. 
Another  disadvantage  is  that  the  reconstructed  wavefront,  at  a  particular 
surface ,  ,  in  Figure  2,  is  a  reproduction  of  the  original  incident  v/ave- 

front  at  the  earth's  surface  only  in  the  case  that  that  wavefront  originated 
from  scatterers  at  a  particular  distance,  depending  on  the  location  of  S 
from  the  earth's  surface.  Observe  that  this  is  a  different  situation  from  that 
that  exists  in  the  reconstruction  of  an  image,  in  which  case,  in  a  given  focal 
plane,  an  image  is  produced  corresponding  only  to  objects  at  given  distance. 

The  difference  in  the  two  situations  is  that  in  the  present  case,  it  is  not  an 
image  of  the  scattering  object  that  is  being  produced,  but  rather  a  reproduc¬ 
tion  of  the  as  yet  unfocused  fields  produced  by  scattering  objects. 

In  Figure  2,  there  is  one  modulator  in  the  modulator  array  for  each 
point  at  which  the  original  field  was  sanqiled.  Beams  of  coherent  light  of 
equal  amplitude  and  phase  are  incident  upon  the  inputs  of  the  modulators. 

Each  modulator  controls  the  amplitude  and  phase  of  the  li^t  incident  upon 
it  in  accordance  with  the  measured  values  of  these  quantities  at  the  receiving 
station  to  which  the  modulator  corresponds.  The  modulator  outputs  pass 
through  the  lens  system,  ,  and  produce  on  the  surface  an 

optical  wavelength  approximation  to  the  microwave  field  that  was  originally 
incident  upon  the  earth's  surface.  This  field  can  then  be  viewed  by  eye, 
just  as  the  virtual  image  produced  by  reconstruction  from  a  hologram  may  be 
viewed.  Alternatively,  if  it  is  desired,  an  optical  system  could  produce  from 
tte  reconstructed  field  a  real  image,  which  could  be  recorded  on  a  photographic 
plate.  The  spatial  locations  of  the  beams  exiting  from  the  modulators  and  the 
phases  of  the  beams  must  be  controlled  to  an  accuracy  corresponding  to  a 
small  fraction  of  a  wavelength  of  the  optical  reconstruction  radiation. 


.1 


It  follows  from  a  spatial  form  of  the  sin  x/x  sampling  theorem  that  the 
interval  between  stations  sampling  the  field  at  the  earth's  surface  must  be 
less  than  or  equal  to  |  X/2L-  ,  which  is  the  same  result  as  was  obtained 


by  a  somewhat  different  method  in  the  preceding  section.  The  relation  be¬ 
tween  array  aperture  size  and  resolution  remains  the  same  as  in  Appendix  1 


Because  it  sqipears  that  the  imesent  method  of  wavefront  reconstruction 
is  more  demanding  (in  that  it  requires  controlled  variation  of  the  phases  of 
the  reconstructing  beams)  than  is  the  reconstruction  scheme  of  Method  1  (or 
Method  2,  for  that  matter),  the  direct  wavefront  reconstruction  method  will 
not  be  discussed  ftirther  here. 
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APPENDIX  3 


SIGNAL-TONOISE  RATIO  ANALYSIS 


In  considering  the  feasibility  of  performing  holographic  imaging  at 
microwave  frequencies  using  the  outputs  of  a  large  coherent  array  of  radar 


sensors,  a  critical  problem  area  is  that  of  the  influence  of  thermal  noise  in 


the  individual  sensors  on  the  signal-to-noise  ratio  (SNR)  of  the  reconstructed 
image.  This  ratio,  for  given  receiver  noise  characteristics  and  a  given  array, 
will  depend  on  the  nature  of  the  object  being  imaged.  This  dependence  is  com¬ 
plex  and  the  effect  of  the  scattering  object  on  SNR  cannot  be  described  solely 
in  terms  of  a  backscattering  cross  section.  It  is  not  possible  to  specify  the 
SNR  that  will  obtain  in  an  image  produced  by  a  coherent  array  by  any  single 
number;  instead,  the  performance  of  the  array  must  be  determined  for  a  range 
of  target  characteristics  sufficiently  broad  to  include  those  of  the  majority  of 
targets  of  interest.  As  a  first  study  of  the  SNRs  resulting  when  holographic 


imaging  is  carried  out  with  such  a  large  aperture  radar  array,  approximate 


expressions  are  developed  in  the  present  analysis  for  the  SNRs  obtained  when 
the  array  is  used  to  image  a  point  (i.e. ,  unresolved)  target,  and  a  target  con¬ 
sisting  of  a  flat  conducting  plate  that  is  assumed  to  be  resolved.  The  method 
of  analysis  is  readily  applicable  to  other  targets.  Results  arc  expressed  in 
terms  of  various  system  parameters,  among  which  are  radar  wavelength, 
transmitter  power,  transmitting  and  receiving  antenna  gains,  coherent  integra¬ 
tion  time,  array  dimensions,  number  of  stations,  and  receiver  noise  tempera¬ 
ture  (assumed  common  to  all  receivers. ) 


PROBLEM  STATEMENT 


For  the  sake  of  definiteness,  a  system  will  be  assumed  consisting 


of  a  single  transmitting  station  illuminating  a  target,  and  a  square  array 
of  dimensions  2L  on  a  side  of  receivers  separated  from  one  another  by 


distance  2b.  It  is  assumed  that  the  real  and  virtual  images  are  sepa¬ 
rated  in  space  from  one  another  and  from  the  other  two  interference 
terms  arising  in  the  holographic  process.  Also,  the  process  described 
in  Ai^ndix  2  as  Method  1  for  reconstructing  a  continuous  image  from  a 
sampled  hologram  will  be  assumed. 

In  order  to  obtain,  with  as  little  delay  as  possible,  an  estimate  of  the 
seriousness  of  the  SNR  problem,  and  of  the  requirements  that  must  be  placed 
on  system  parameters  in  order  to  obtain  an  acceptable  SNR,  a  number  of  ap¬ 
proximations  have  been  used.  These  will  be  pointed  out  as  they  are  intro¬ 
duced.  In  the  future,  a  more  exact  analysis  of  the  problems  considered  here 
will  be  conducted,  Insofar  as  such  treatment  seems  required. 

NOISE  ANALYSIS 

Coherent  Integration- Equivalence  of  Receiver  Noise  to  an  Incident  Field 

It  is  assumed  that  each  receiving  station  is  arranged  to  receive  a 
scattered  field  of  the  same  polarization.  If  E^  is  the  complex  field  intensity 
(of  the  polarization  component  of  interest)  incident  on  a  receiving  antenna, 
the  power  flux  per  unit  area  carried  by  the  wave  is  1 E^  |  ^  /2t)  ,  where  tj  is 
the  impedance  of  free  space.  If  the  gain  of  the  receiving  antenna  is  , 
the  power  absorbed  by  a  matched  load  attached  to  the  antenna  is 


P,=  \ES 


(1 
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If  the  antenna  and  load  are  repreuented  by  the  circuit  of  Figure  1, 


where  R  is  the  radiation  resistance  of  the  antenna,  P  may  be  written 


(2) 


where  V  is  understood  to  be  a  peak  value.  Thus 
a 


and 


\\/^\  »  aVz  P»  Ra  , 


(2) 
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If  a  sequence  of  N  pulses  is  received  from  a  given  target  and  added  coherently, 

the  magnitude  of  the  resultant  voltage  will  be  N  |  V  I  ,  if  I V  |  is  that  corres- 

a  a 

ponding  to  a  single  pulse.  Further,  by  Equation  (4)  the  field  intensity  that 
would  lead  to  the  voltage  n]v  |  is  Just 

^IL  (5) 


Now  consider  the  noise.  Assiune  that  the  background  noise  picked  up 

from  the  sky  by  the  receiver  is  characterized  by  a  temperature  T  ,  the 

B 

noise  due  to  antenna  losses  by  temperature  ,  and  the  noise  introduced  by 

the  input  stages  of  the  receiver  by  tenoperature  T  .  The  total  noise  tempera- 

xl 

ture  is  then  T  =  T  +  T.  +  T  .  If  the  bandwidth  of  the  receiver,  Af  ,  is 

xS  A  li 

assumed  to  be  small  compared  with  the  center  frequency,  f  ,  the  resulting 
noise  may  be  pictured  as  a  carrier  wave  of  frequency  f  of  amplitude 
and  phase  ^  ,  where  A^^  and  are  functions  of  time  that  cannot  vary 
substantially  in  time  intervals  of  duration  less  than  1/Af.  Assuming  that  the 
receiver  bandwidth  is  approximately  that  corresponding  to  a  received  pulse  of 
duration  Tp  (i.  e. ,  Af  ~  1/Tp)  ,  the  effect  of  the  noise  is  to  add  a  voltage 
vectorially  to  the  vector  voltage  of  each  single  pulse.  Thus,  if  for  the  jth 


pulse  the  received  voltage  is  v  ,  ,  and  if  the  complex  noise  voltage  during 

rj 

the  interval  of  that  pulse  is  n^  (where  n^  can  be  considered  constant  over 
the  duration  of  a  sirgle  pulse)  ,  the  total  voltage  recorded  by  the  receiver  for 


the  jth  pulse  is 
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"^1  =  x-j  +  ''tj 


where  the  magnitude  of  is  that  of  the  (substantially  constant)  envelope  of 
the  noise  during  the  time  of  the  jth  pulse,  and  the  phase  of  m  is  the  (sub¬ 
stantially  eonstant)  phase  of  the  noise  during  that  pulse.  The  sum  of  N  sueh 
pulses  is 


=  2  4  Mj')  =  N  yir  2  ">^1 


where  it  is  assumed  that  v  .  is  eonstant  with  value  v  (that  is,  it  is 

rj  r 

assumed  that  the  sueeessive  reeeived  pulses  are  phase  shifted  so  that  they 
may  be  added  in  phase).  Then, 


iXil*  «  (Ny;.  +  2'ns)(N -K?  +  E '<) 


=  N*"  \  y:rV‘  +  N 


where  |ni  =  ^ J*  The  bar  indieatcs 

an  ensemble  average. 
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According  to  Equation  (4),  a  voltage  at  the  input  of  the  receiver  of  mean* 
square  magnitude  N  |n|  would  be  produced  by  a  single  incident  pulse  having 
electric  field  of  mean  square  magnitude 


1"  =  UL 

However,  the  mean  square  noise  amplitude,  for  the  case  of  narrow  band  noise 
is  (this  is  twice  the  mean  square  noise  voltage**), 


|/nl^  =  S  fe  T 


(7) 


SO 


(8) 


The  noise  amplitude  and  phase  to  be  added  to  a  given  pulse  at  a  particu¬ 
lar  station  is  considered  to  be  independent  of  those  quantities  for  all  pulses, 
including  the  corresponding  pulse,  at  all  other  stations  (this  may  not  be  true 
for  the  noise  due  to  the  sl^  background),  and  to  be  independent  of  those 

♦The  word  "mean"  will  be  used  consistently  to  refer  to  an  ensemble  average. 

**The  noise  voltage  is  assumed  to  be  Gaussian  with  zero  mean.  Under  this 
assumption,  the  envelope  of  narrow  band  noise  is  governed  by  a  Rayleigh 
distribution,L^^l  and  it  is  strai^tforward  to  show  that  the  mean 
square  value  of  this  envelope  is  twice  the  mean  square  instantaneous 
noise  voltage. 
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quantities  for  all  other  pulses  at  the  given  station.  The  pulse  amplitudes  will 
be  Rayleigh  distributed  for  the  narrow  band  noise  case,  and  the  phases  dis¬ 
tributed  uniformly  over  (0,  2?:)  .  It  is  assumed  that  at  each  station  the  inter¬ 
val  between  pulses  is  long  compared  to  1/Af.  Thus,  the  noise  voltage  that 
is  present  at  each  station  after  the  coherent  integration  of  a  train  of  N  pulses 
may  be  represented  by 


and  is  Rayleigh  distributed,  llie  quantities  and  (p  are  independent 
N  N 

at  the  different  stations,  and  are  independent  of  each  other  at  a  given  station. 
(As  indicated  above,  Ej^  and  tp  at  different  stations  may  not  be  independent, 
it  the  predominant  noise  component  arises  from  the  sky  background. ) 

Noise  Terms  in  Hologram 

If  he  reference  signal  used  in  forming  the  hologram  is  E  (x)  ,  where  x 

K 

is  a  two-dimensional  vector  describing  a  point  on  the  hologram  surface,  the 
total  field  falling  on  the  hologram  plane  is 

Et.ix)=  E6<5'^  +  E.^cx)-t  (10) 


Here  E_(r.)  is  the  field  scattered  by  the  object  to  be  imaged,  and  E  (x)  is 
D  N 

the  fictitious  noise  field  described  in  the  preceding  section.  The  power  flux 
falling  on  the  hologram  plane  is  then  \e^\^ / 2 r}  .  The  squared  term,  |e.j,1  , 
may  be  written  in  full  as 
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IEtI  “  lEol*"  +  lEnl*' +  lEiil’’ +  £■>£«*  t  ErEo 

+  e,e;  +  e,e: +e,e*^EheS. 

If  the  signal,  E,j,  ,  is  used  to  form  a  hologram  (either  continuous  or  sampled), 
and  an  image  is  reconstructed  from  the  hologram  by  Method  1  of  Appendix  2,  the 

f  .  and 

process  that  are  described  in  Appendices  1  and  2.  In  particular,  the  terms 

2 

and  lead  to  the  production  of  images,  and  the  terms  |Ej^| 

and  I E  I  ^  to  undesired  components.  By  suitable  placement  of  the  reference 

IV 

signal,  it  can  be  arranged  that  the  fields  produced  in  the  reconstruction  by  the 

terms  |e  |  ,  |e  |^  ,  and  EE*  ,  are  not  brought  to  a  focus  in  the  same 
U  K  K  U 

plane  as  the  field  due  to  E  E  *  .  Experience  gained  in  optical,  on- axis, 

D  K 

Fraunhofer  holography  indicates  that  weak,  out-of-focus  fields  do  not  seriously 
degrade  the  desired  image.  The  effects  of  the  terms, 

IE  E,E„%E„E^E,E:tE„E;f^  (U) 


Ej^Ejj  lead  to  results  in  the  reconstruction 


terms  |Ej^|  ,  |Ej^ 


involving  the  thermal  noise,  will  now  be  investigated. 

Effect  of  Noise  Terms  in  Reconstruction 

According  to  Equation  (6)  in  Appendix  2,  the  noise  terms  of  Equa¬ 
tion  (12)  will  give  rise  to  the  following  sum,  S ,  which  must  be  evaluated 
in  order  to  determine  the  effect  of  these  terms  on  the  reconstructed  image 
of  the  target  object. 
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•  -L^  ('»n«j,m) ,  (13) 


The  various  symbols  appearing  in  Equation  (13)  have  been  defined  in  the 
preceding  Appendices.  The  mean  squared  magnitude  of  S  is 

ji/*t  - 

Note  that  I..(x)  is  real.  This  has  been  used  in  Equation  (14).  In  order  to 
N 

evaluate  S  S*  ,  write 


Here  the  subscripts  j,  n  indicate  the  evaluation  of  a  quantity  at  x  -  mu. 

J*n 

I.et  a.  =  E„.  +  E_.  .  Then  Equation  (15)  becomes 

j,n  Dj,n  Rj,n 


Now  observe  that 
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because  the  amplitude  and  the  phase  of  ^  are  independent,  and  the  phase 

is  uniformly  distributed  over  (0,2ir),  E„,  =  0.  Using  this,  and  observing 

NJ,  n 

tl^t  a.  is  not  random,  it  is  seen  that 


Substituting  Equations  (16)  and  (18)  into  Equation  (17)  and  rearranging,  there 
results 


Here  the  noise  E^^  ^  is  assumed  statistically  homogeneous,  so  that,  for 

example,  |Ej^j  1^  .  The  independence  of  the  noise  at  the 

different  receiving  stations,  i.  e. ,  for  different  pairs  J,  n,  is  now  used  to 
simplify  Equation  (19) 

“  C  ■  lEhi-i  +  Ehj,- 

+ i  TEZx  *  4.1!^  h,,.,. 
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However,  again  using  the  independence  of  amplitude  and  phase  of  E„.  and 

Nj,n 

the  uniformity  of  phase  distribution  over  (0, 2ir )  , 


and 


E  ’  =  e;.  -o 


Hence 


+  i^K  r 


(20) 


Equation  (20)  can  be  simplified  somewhat  further  by  using  the  fact  that 

!  E^,  1  is  governed  by  a  Rayleigh  distribution  to  write 

•  Nj,n 


Using  statistical  homogeneity  to  drop  the  subscripts  j,n  from  E^^  , 
Equation  (20)  becomes 


(21) 
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Now  substitute  Eq,iaUM*(2lT^St?1E5Suior^^^'^[n^^umo^j^n^o 


obtain 


.  ss*»  2  lE«r  (IEhI'+iK^i*') 


JiMl 


=  N  I  EmP  i-  2  IE«r  +  KV-  P 


(22) 


where  M  is  the  total  number  of  stations,  and 

P=  12. 

I  1 


I*- 

Ca.  I 


(23) 


In  order  that  the  reconstructed  fields  corresponding  to  the  terms  |  E  ] 

,  ,2  ♦  ^ 

|Ej^  I  ,  aud  ,  will  not  be  in  focus  in  the  same  plane  as  that  in  which 

the  field  due  to  E  E_  is  in  focus,  it  is  desirable  that  the  phase 

li 

I  12 

k'g  |u  I  /2  vary  over  the  aperture  (i.  e. ,  over  j,  n )  by  an  amount  at 

J>n 

least  equal  to  2v  ,  and  perhaps  by  considerably  more.  Considering  also  the 

phase  change  over  the  range  of  J,n  due  to  the  term  k'  a  •  u  /  ^9  > 

j,  n  L 

appears  that  the  magnitude  of  the  siun  in  Equation  (23)  can  be  estimated  by 
considering  the  sum  to  represent  a  random  walk  in  the  complex  plane,  each 
step  being  of  unit  length.  Thus,  as  a  crude  estimate,  P  ~  M  .  This  esti¬ 
mate  will  be  refined  in  future  work.  For  the  present,  however,  estimate 
P  by  M  and  substitute  in  Equation  (22)  to  obtain 
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2  1eJ^(m  I£J'  -h  'Z\«.i..\') 


Equation  (24)  Involves,  in  the  quantity  a,  ,  the  strengths  and  the  dis- 

j.n 

tributions  over  the  hologram  plane  of  the  field  scattered  from  the  target 
and  the  reference  field  E^  .  Using  the  definition  of  > 


S  Ui.-r  *  2  )Ed.  +E..  1 


A  careful  evaluation  of  this  summation  would  be  tedious;  however,  for  present 
purposes,  it  is  probably  sufficiently  accurate  to  assume  that  the  contribution 
to  the  summation  of  tlte  middle  two  terms  in  the  above  egression  is  small, 
because  of  the  variation  in  phase  of  those  terms  in  the  summation  over  J 
and  n  .  If  this  assumption  is  made,  the  above  expression  becomes 


I." 


in  which 


Eo”;  - 


-Jr  2  lEti-r 
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and 


(27) 


2  2  2  2 
Thus,  E  .  and  E  are  the  average  values  of  |e_,  |  and  ]e„,  |  , 

DA  RA  '  Dj,n  Rjin' 

respectively,  over  the  various  receiving  stations.  These  averages,  in  turn, 

are  approximately  equal  to  averages  taken  over  the  entire  area  of  the  hologram. 

Using  Equaticm  (25),  Equation  (24)  becomes 


SSSaM  +£.*). 


(28) 


2  2 

Consider  now  the  averages  E  .  and  E  .  .  Two  cases  will  be 

1a\  RA 

studied:  a  point  target  (i.e. ,  an  unresolved  target),  and  a  flat  plate  (assumed  tu 
be  resolved). 


Case  i  ~  Point  Target 

From  the  first  term  of  Equation  (7)  of  Appendix  1,  the  signal  field  from 
a  point  source  located  at  is  [taking  D(i)  A*  6(1  -  (g)  ,  with  A' 
constant]  , 


Here  the  quantity  x^  ,  appearing  in  Appendix  1,  is  taken  to  be  zero.  The 
c(mstant  A'  will  be  evaluated  by  equating  the  energy  flow  at  x  =  0  predicted 
by  Equation  (29)  to  that  that  would  exist  at  the  location  if  there  were  a  point 
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radar  reflector  at  |  scattering  energy  from  some  illuminating  signal  toward 
the  point  X  =  0  in  the  hologram  plane.  L  ^  a  be  the  scattering  cross-sec¬ 
tion  describing  this  process.  Then 

•l.  A  »  ‘ 


kM'  _  Pr<ST<r 


and 


i/T 


(30) 


where  is  the  transmitter  power  (during  each  individual  pulse),  is 
the  transmitter  antenna  gain,  and  f  ^  is  the  range  from  the  transmitter  to 
the  scattercr  (£^  «  £).  Substituting  Equation  (30)  into  Equation  (29),  and 
multiplying  by  N  to  account  for  coherent  integration  of  N  pulses. 


!i 


,<t (i * !2l!  -  4  l^l*) 

m  IV  ^  ■  _  1'  ' 


E.c5i)=- e'^'  T~  ■  iT'  M 


The  reference  signal  is  obtained  from  the  second  term  of  Equation  (7) 
of  Appendix  1  by  setting  6  =  X3  =  0  .  Thus, 


E,(7) 


B  ik[^s«+ 


1. 


(32) 


3R 


in  which  B  is  a  constant,  and  the  other  terms  are  defined  in  Appendix  1 , 
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To  evaluate  Ej^  and  Ej^  of  Equations  (26)  and  (27),  inasmuch  as  u 

point  scatterer  cannot  give  rise  to  fine  structure  in  its  scattering  pattern,  it 

12  2 
E^.  1  and  |e^  I  are  con- 

Dj,n'  '  Rj,n' 

stant  over  the  receiving  array.  In  this  case,  using  Equations  (31)  and  (32), 


‘hA 


a.  ; 


(33) 


and 


E" 


|Bl^ 


a  +  T,,)" 


(3-1) 


Now,  turning  to  Equation  (6)  of  Appendix  2,  observing  th(3  definition  of 
the  sum  S  in  Equation  (13),  it  is  apparent  that  the  mean  squared  magnitude 
of  the  contribution  to  the  reconstructed  field  of  the  noise  terms  in  the  holo¬ 
gram  may  be  written 


IE  ss*=k:s5 

n 


(33) 


where  C  ,  f  ,  z  ,  z  ,  andm  ,  are  constants. 

X  M 

In  obtaining  Equation  (35)  from  Equation  (6)  of  Appendix  2,  the  sinx/x 
factors  in  Equation  (6)  have  been  ignored,  because  any  constant  multiplication 
factor  arising  from  these  factors  would  also  appear  in  the  expression  that  will 
be  obtained  in  the  following  work  for  tliC  squared  field  intensity  of  the  image  of 
the  target.  Thus,  when  the  ratio  of  the  squared  signal  field  to  the  mean 
squared  noise  field  is  taken,  the  factors  will  cancel.  It  should  also  be  observed 
that,  although  Equation  (6)  is  written  with  particular  reference  to 
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reconstruction  of  an  image  from  the  hologram  term  I  ,  the  expression  of 

2 

Equation  (6)  is  equally  sq^plicable  to  the  calculation  of  the  field  appearing  in 
the  reconstruction  arising  from  any  other  term  in  the  hologram. 


Case  ii  -  Flat  Plat  Target 


The  first  term  of  Equation  (7)  of  Appendix  1  gives  for  the  field  at  the 
hologram  surface  due  to  the  target 


-  jii 


D(f)e 


!ll^) 

) 


(36) 


where,  as  before,  x  has  been  taken  as  zero.  Assume  the  disturbance  D(4) 
giving  rise  to  the  field  to  be  of  the  form 


in  which  B'  is  a  constant  (assumed  real)  and,  if  \  ~  »  where 

XX  6 

1  and  i  are  orthogonal  unit  vectors  parallel  to  the  image  plane, 

X  M 


I  ^  ^  i  Yi  1  — eiiTjici 


Assuming  the  hologram  surface  to  be  in  the  Fraunhofer  region,  the  term 
k\l\^  /2i  in  the  exponent  in  Equation  (36)  will  be  neglected.  Substituting 
for  D(i)  in  Equation  (36)  and  integrating  yields, 
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in  which  x  =  l^x^  +  constant  B'  may  be  expressed  in  terms  of 

radar  parameters  and  parameters  describing  a  familiar  target  by  observing 
that,  if  the  radar  waveleng^  is  much  less  than  2a,  the  backscattering  cross 
section  of  a  perfectly  conducting,  square,  flat  plate  of  sides  2a  is 


>s^ 


(38) 


Hence  at  x  -  vO,  equating  the  power  flow  that  Equation  (37)  predicts  will 
strike  the  hologr  am  plane  to  that  that  would  strike  the  plane  in  the  case  of 
scattering  from  a  flat  plate  with  backscattering  cross  section  cr^  , 


^  Pr  Gt 


/ 


Pt  (^T 
27r  Jlj 


(39) 


Substituting  Equation  (39)  into  Equation  (37)  and  multiplying  by  N  to  account 
for  coherent  integration  of  N  pulses,  there  is  obtained 


K 


t 


Sim. 


(40) 
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The  average  ^  can  be  approximated  by  a  continuous  average  of  1  E^(x)l^ 
U/V  u 

over  the  region  of  the  hologram,  i.  e. ,  over  -  L  ^  ^  L  and 

-L  ^  X  ^  L,  The  result  is 
2 


E 


BA 


N'  Pt 


(41) 


2 

The  expression  for  E  is  given  by  Equation  (34).  Equation  (35)  again 

It  A 

gives  the  mean  squared  value  of  the  noise  field  in  the  reconstructed  image, 

2 

provided  that  the  expression  for  E^,  given  in  Equation  (41)  is  used  in 
evaluating  SS  in  Equation  (35). 

INTENSITY  OF  RECONSTRUCTED  SIGNAL** 

A  point  target  and  a  flat  plate  target  are  again  considered  separately. 

Point  Target 

In  this  case,  in  the  reconstructed  image  the  target  is  not  resolved. 
There  will  be  the  diffraction  pattern  of  the  hologram  aperture  where  the 
point  image  should  appear.  The  squared  field  value  at  the  center  of  this 
pattern  will  be  evaluated  and  compared  with  the  mean  squared  noise  field, 
in  order  to  determine  the  SNR.  The  actual  diffraction  pattern  resulting  by 
reconstruction  from  a  sampled  hologram  will  differ  somewhat  from  that 
obtained  by  reconstruction  from  a  continuous  hologram,  but  for  present  pur¬ 
poses  this  difference  will  be  neglected.  *  In  the  case  of  a  continuous 

*  The  result  obtained  by  this  assumption  (i.  e. ,  Equation  (46))  is  the  same  as 
that  that  would  be  obtained  by  considering  in  detail  the  reconstruction  from  a 
sampled  hologram,  provided  that  the  number  of  receiving  stations  is  large 
compared  to  unity. 

See  addendum,  p.  161. 
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hologram,  from  Appendix  1, 

r  .  -  /i^  k  ^  u  ) 

In  obtaining  Equation  (42)  from  Appendix  1,  the  displacement  vector 
has  be^^n  set  equal  to  zero,  and,  for  the  fimction  D  ,  the  constant  value 
(2ir)  ^  has  been  substituted  in  accordance  with  Equation  (19)  of  Appendix  1 
(1.  e. ,  a  point  scatterer  has  been  assumed).  The  quantity  F  is  a  constant. 
If  the  original  hologram  wa«  square,  of  length  2L  on  a  side,  and  if  the 
hologram  was  scaled  down  by  a  factor  m  before  optical  reconstruction,  the 
above  integral  should  be  taken  over  a  square  region  centered  at  the  origin 
and  length  2L/m  on  a  side. 

The  result  of  the  integration  is*' 


L  ( 
-  I 


J*  hn 


Tja  / 


.\ 


Thus,  in  the  case  in  which  the  source  term  D(|)  in  Eqxiation  (36)  is  of  the 
form  6{l)  ,  in  the  reconstruction  (he  squared  field  at  the  point  at  which  the 
image  of  the  target  should  appear  is 
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SIGNAL  TO  NOISE  RATIO 
General  Result 

Using  the  foregoing  results,  it  is  possible  at  once  to  obtain  the  ratio 
of  the  squared  magnitude  of  the  signal  component  of  the  image  field  to  the 
mean  squared  magnitude  of  the  noise  component. 


Dividing  Equation  (46)  by  Equation  (35),  and  substituting  from  Equatiunia 
(8),  (28),  (33),  and  (34),  there  is  obtained,  after  some  simplification, 


l£r«t 


Pr  CrrCn 

31/r  +  /■,<)' 


r^r./vTTA-f  /  Pr  CT<r>/  ^ 

*L 

(4ri) 

Some  implications  of  Equation  (47),  and  some  simplification  of  that 
equation  to  special  cases  will  be  considered  later.  First,  however,  the 
flat  plate  target  will  be  treated. 

Flat  Plate  Target 

Divide  Equation  (47)  by  Equation  (35),  and  substitute  from  Equations 
(8),  (28),  (34),  and  (41)  to  obtain,  after  simplifying 
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.rr;  l*  ‘ ^ 
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simplification  and  Numerical  Results 


For  brevity,  let 


I 

i  and 


M«ii< 


(49) 


/e^,. 

'll*  it 

2  2  I  i2 

Equations  (48)  and  (49)  indicate  that  and  start  at  zero  for  |  B| 
equal  to  zero,  and  increase  monotonically,  asymptotically  approaching  con¬ 
stant  values,  as  ]  b]  ^  becomes  larger.  The  value  of  |  b]  ^  that  is  required 
2  2 

to  bring  R  or  R  to  a  value  that  is  one-half  of  its  asymptotic  maximum 

O  u 

is  that  value  of  |b1^  for  which  the  quantity  2lBlV(f  +  is  equal  to 

the  sum  of  the  other  terms  in  the  factor  in  brackets  in  the  e3q)ression  in 
question,  i.  e. ,  a  value  of  reference  signal  power  roughly  equal  to  the  sum  of 

1  i2 

the  power  of  the  target  signal  and  the  noise.  Increasing  (  B|  beyond  this 
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valttt  brings  R  ^  or  R-^  still  closer  to  the  maximum,  with  a  value  of  R  ^ 

2  ^  2  ® 
or  R^  of  0.  d  of  the  maximum  ooourrlng  when  |  b|  is  nine  times  the 

abovementlcned  value. 

Ibe  asymptotlo  maximum  SNR  is  found  from  F^iuation  (48)  or  Equation 
2 

(49)  fay  letting  ]  b|  iq>proach  infinity.  Hie  result  is: 


L  //  f  r  ^  4k 


1*  ^  ^  ^ 


H  ’  _ _ '  '  /I  »»  ir  wt  ^ 


v'-A^rrtrC,H 

ff  tile  bandwldtii  of  the  receiver  is  matched  to  the  pulse  bandwidth,  the 
poise  duration,  Tp  ,  and  the  bandwidth,  A  f ,  a^iiearinf  in  Equations  (50) 
and  (51),  are  related  by  ~  i/Af  . 

Daflning  tiie  dimensionless  "radar  parameter" 


and 


(53) 


3V>(;-b‘'A1  (64) 

Note  alao  that,  by  Equations  (46)  and  (47),  the  magnitude  of  the 
reconstructed  field  at  the  center  of  the  dlffractimi  spot  for  a  point  target 
will  equal  the  magnitude  of  the  reconstructed  field  in  the  interior  of  the 
image  of  a  plate,  if  the  scattering  cross  section  of  the  point  target  is 


HT- 

(55) 

2  2 

This  value  of  of  a  also  corresponds  to  equality  of  R  and  R, 

^  o  max  L  max 

As  a  numerical  example  ,  consider  the  case  in  which  the  various 
parameters  have  the  following  values: 

P_  ■  10®  watts, 

T 

>  5  X  10®  (~57db)  ~  90  foot  dish  at  S-band, 

N  -  1, 


r 

i 
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and 


X  =  0.  Im  (S-band) 

=  I  =  3  X  lO^m, 

T  =  300°K, 

Tp  =  l/Af  «  10  ®  sec, 

4 

L  s  1.5  X  10  m, 

M  =  121, 

-2  2 

a  =  10  m  , 

b  =  1.50  X  lO^m, 

-»23  ^ 

k  =  Boltzmann's  constant  =  1. 38  x  10  joule/°K. 
The  radar  parameter  is 


24 

F_  -  7. 82  X  10  , 

it 

R  ^  =  0. 32  or  -4. 9  db  , 

o  max 


R,  max  »  4. 0  X  10  or  46  db 


ahd 

2 

»  1260m  or  31  dbsm. 

The  above  SNR  for  the  plate  Is  quite  large,  but  the  ratio  is  low  for 
the  point  source.  If  it  were  desired  to  increase  the  SNR  1^,  say,  20db,  so 
that  point  sources  of  low  cross  section  would  be  more  clearly  visible,  it 
would  probably  be  necessary  to  use  coherent  integration.  Thus,  in  the  above 
case,  if  N  had  been  100,  the  result  would  have  been 


158 


M6i)-2 


N  =  100 


o  max  =  15  db 


2  max  =  66  db  . 

The  ^eat  difference  in  SNR  for  the  point  source  and  the  plate  may  indicate 
that  some  form  of  amplitude  compression  would  be  advisable. 

If  new  holograms  were  required  of  the  target  complex  at  a  rate,  fj^  , 
of  100  per  second,  for  example,  the  average  power  required  of  the  trnnsmitter 
would  be 


P  =  P.„T„Nf,_ 
a  T  P  b 


=  10  watts,  for  N  =  100. 


Regarding  the  coherent  integration  time,  ,  note  that  this  interval 

must  be  sufficiently  short  so  that  in  the  time,  ,  no  two  parts  of  the  target 

complex  that  have  significant  cross-section  can  move  more  than,  for  example, 

a  tenth  of  the  operating  wavelength  relative  to  one  another.  If  10  meters/sec 

is  taken  as  a  reasonable  upper  limit  to  the  relative  velocities  of  such  sc.at- 

-3 

tcrers,  the  maximum  permissible  value  for  T  is  10  sec.  If  the  pulse 
•6  ^ 

duration  is  10  sec,  and  N  =  100,  it  is  seen  that  the  total  transmitting 

-4 

time  required  for  producing  a  single  hologram  is  10  sec.  The  100  pulses 

required  for  a  single  hologram  may  be  easily  fitted  within  the  maximum 

-3 

permissible  time  interval  of  10  sec.  It  is  understood  that  in  performing 
the  coherent  integration,  compensation  is  made  for  the  rather  uniform 
phase  change  that  is  occurring  from  pulse  to  pulse  due  to  the  translational 
motion  of  the  target  complex.  The  use  of  an  individual  pulse  duration  of 
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10  sec  prevents  significant  phase  change  from  occurring  due  to  target 
translational  motion  during  a  single  pulse.  Of  course,  rather  than  a 
sequence  of  discrete  one-microsecond  pulses,  a  single  long  pulse  could  be 
used,  making  appropriate  compensation  for  phase  change  due  to  translation. 
(The  effect  on  the  results  of  the  present  analysis  of  the  use  of  a  single  long 
pulse  will  be  investigated.  Some  changes  in  the  results  may  occur,  because 
it  has  been  assumed  that  the  noise  on  successive  pulses  is  independent. ) 

In  the  numerical  example  presented  above,  the  requirements  on  the 
.  receiver  stations,  aside  from  that  of  coherent  integration,  have  been  made 
fairly  modest;  i.  e. ,  an  antenna  gain  of  10  and  a  noise  temperature  of  300^K. 
lliat  the  individual  receiving  station  be  as  simple  as  possible  is  dictated  by 
the  required  number  of  such  stations  (121  in  the  example).  On  the  other 
hand,  only  one  (or  at  most  a  few)  transmitter  is  needed,  so  that  the  require¬ 
ment  of  a  90  fcot  dish  at  the  transmitter  is  not  imreasonable. 

In  the  future,  continuing  attention  will  be  given  to  the  problem  of  SNK. 
In  particular,  the  effect  of  reccmstructing  the  sampled  hologram  by  methods 
other  than  the  one  that  has  been  assumed  here  will  be  investigated. 
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ADDENDUM 

On  pages  151  and  152  of  M68-  2,  in  obtaining  the  intensity  at  the  center 
of  the  diffraction  spot  resulting  from  the  reconstruction  of  the  image  of  a 
point  target  of  unit  strength  (i.  o. ,  a  6-fxmction, )  there  was  used  for  simplicity 
the  expression  tliat  had  been  obtained  for  reconstruction  of  a  continuous  hologram. 
In  the  event  that  b'  =  b  (i.e. ,  the  holes  in  the  scaled  down  sampled  hologram 
arc  so  large  that  they  abut, )  such  a  procedure  is  legitimate:  however,  when 
b’  f  b  ,  it  is  necessary  to  use  the  reconstruction  expression  for  a  sampled 
hologram.  If  this  is  done,  using  Equations  (6)  and  (7)  of  Appendix  2,  and 
assuming  the  focusing  condition  C  , 

4t 


4F  L 

E^^  («') - 

rec  ^ 

m 


.,2  'b*  k'a  'b' 

b'  1  2 

sine  -  sine  - 

”“2  ”"2 


where 

"^2<R 

the  location  of  the  center  of  the  diffraction  spot.  The  above  expression  for 
E  (a*)  replaces  Equation  (43)  of  Appendix  3  of  M68-2.  Neglecting  sine 

170C 

function,  as  explained  on  page  148.  Equation  (44)  of  Appendix  3  is  replaced 

2  4 

,2  b''* 


Equations  (45)  and  (46)  of  that  Appendix,  were,  in  fact,  derived  using  the 

correct  expression  for  E  (ot*)  .  They  remain  unchanged. 

rec 
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